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20.  ABSTRACT 

This  report  is  in  two  parts.  In  the  first  part  averaging 
projections  on  finite  dimensional  subspaces  are  used  in  order 
to  approximate  an  autonomous  nonlinear  functional  differential 
equation  with  state  space  RnxLP  [-r,0;Rn)  by  a sequence  of 
ordinary  differential  equations.  This  type  of  approximations 
already  proved  to  be  successful  for  the  numerical  treatment  of 
hereditary  control  problems  involving  linear  functional 
differential  equations. 

In  the  second  part  the  authors  investigate  the  optimal  control 
problem  for  Volterra  integrodif ferential  equations. 

X(T)  = A(t)x(t)  + ^ F (t-sjx(s)ds  + B(t)u(t) 

where  the  target  sets  are  elements  of  some  function  space, 

X(t)  =^(t)  for  0<T-r£t£T 


The  approach  used  is  the  abstract  theory  of  Dubovitskii- 
Milyutin,  the  result  is  a necessary  condition  in  form  of 
a maximum  principal. 


Projection  Methods  .in  Nonlinear  Autononuous  Functional 
Differential  Equations. 


Abstract 

In  this  part  averaging  projections  on  finite  dimensional  sub- 
spaces are  used  in  order  to  approximate  an  autonomous  nonlinear 
functional  differential  equation  with  state  space  RnxLp  [ - r , 0 ;H!n  ) 
by  a sequence  of  ordinary  differential  equations.  This  type  of 
approximations  already  proved  to  be  successful  for  the  numerical 
treatment  of  hereditary  control  problems  involving  linear 
functional  differential  equations. 
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& 1.  Introduction 

For  the  numerical  treatment  of  hereditary  control  problems 
various  authors  have  develloped  an  approach  which  uses  pro- 
jections on  finite  dimensional  subspaces  of  the  given  state 
space  (see  especially  [4]  and  [2]  and  the  extension  discussion 
on  related  literature  given  in  [2].  Problems  considered  are 
of  the  following  type: 

(P)  Minimize  $(u)  = J(x(-,u),u)  subject  to 

x(t,u)  = f (x(t,u) ,xt(u)  ,u(t) ) , telO.tj], 
x ( 0 , u ) = n£Rn, 

x(s,u)  ^ <p(s)  a.e.  on  [-r,0),  tp€Lp( -r  , 0 ; Rn ) , 

and 

ueu,  U a closed  convex  subset  of  L^(  [0,t  ^ ] , IRm)  . 

For  any  function  x:[-r,a)  -»  Rn  , a > 0,  the  symbol  xt  denotes  the 

map  s -*  x(t  + s),  se[-r,0],  t > 0.  The  state  space  used  in  (P) 

is  X = RnxLp(-r,0;IRn)  and  is  therefore  infinite  dimensional. 

A possible  approach  to  problem  (P)  can  be  described  as  follows. 

N 

Choose  a sequence  (X  ) of  finite  dimensional  subspace  of  X such 

N N+l  N N N 

that  X c X for  all  N with  projection  n : X -*  X , it  -»  idY 

N A 

as  N -»  oo.  Define  operators  G such  that  the  solution  of  the 

functional  differential  equation  in  (P)  is  approximated  by  the 
N 

sequence  (x  (t;u))  of  solutions  of 

iN(t)  = GN(xN(t),u(t))  , t > 0, 
xN(0;u)  = nN(n,<p)  . 

N N N 

Moreover,  we  have  to  define  a sequence  4>  (u)  = J (x  (-.uj.u) 
of  cost  functionals. 

Therefore  in  addition  to  (P)  we  have  a sequence  of  problems: 


(PN)  Minimize  $N(u)  subject  to 
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xN(t;u) 
xN(0;u) 
a nd 
u£U. 

N - - N 

Suppose  that  (P)  and  (P  ) have  the  solutions  u,u  , respectively. 

- N 

Then  the  question  is:  Under  what  condition  we  have  u -»  u (weakly 
or  strongly),  xN(t;uN)  -»  x(t;u)  and  $N(uN)  -►  <t>(u) 

If  convergence  can  be  prooved  then  we  have  an  approximation  of 

the  infinite  dimensional  problem  (P)  by  the  sequence  of  finite 

N 

dimensional  problems  (P  ). 

The  approach  indicated  here  was  applied  successfully  for  the 
numerical  solution  of  problems  (P)  where  the  functional  differen- 
tial equation  is  linear  (cf.  [2],  [3]).  The  following  sections 
deal  with  the  approximation  problem  for  the  functional-differential 
equation  only  assuming  local  Lipschitz  type  conditions.  The  results 
presented  here  constitute  the  first  but  essential  step  towards 
the  numerical  treatment  of  nonlinear  hereditary  control  problems 
by  projection  methods. 

As  mentioned  above  an  important  feature  of  this  investigation 
is  that  only  local  Lipschitz  type  conditions  are  assumed.  This 
implies  that  in  general  we  only  have  local  existence  of  solutions. 
But  even  in  the  case  where  we  have  global  existence  of  solutions, 
the  solution  semigroup  in  general  is  not  a semigroup  with 
dissipative  infinitesimal  generator.  This  is  only  the  case  in 
special  situations  [20).  Therefore  the  results  of  the  existing 
theory  of  nonlinear  semigroups  cannot  be  applied  directly. 

In  Section  2 we  provide  the  existence,  uniqueness  and  continuous 
dependence  results  which  we  need  in  this  paper.  Section  3 gives 
an  analysis  of  the  local  semigroup  generated  by  solutions  of  the 
underlying  f unc t i ona 1 - d i f f eren t i a 1 equation.  The  main  idea 


= GN(xN(t;u) , u ( t ) , t£ [ 0 , t i ] , 

N,  . ' 1 

= n (n  ,if>)  , 


in  Section  4 is  to  approximate  the  solution  semigroup  by  semi- 
groups with  nicer  properties  (i.e.  dissipative  infinitesimal 
generator).  These  approximating  semigroups  are  solution  semi- 
groups of  "averaged"  functional  differential  equations.  In 
Section  5 the  results  of  Section  4 are  carried  over  to  the 
local  case.  Finally,  in  Section  6 the  averaging  projections 
introduced  in  this  connexion  by  Banks  and  Burns  [2]  (see  also 
the  extensive  literature  given  there)  are  considered  in  the  case 
of  difference-differential  equations  where  the  right-hand  side 
is  locally  Lipschitz.  For  Lipschitz  continuous  initial  data  we 
get  uniform  convergence  of  the  approximating  ordinary  differential 
equation  solutionson  compact  intervals. 

The  methods  applied  in  the  following  sections  will  be  carried 
over  to  time-dependent  equations  and  equations  with  infinite 
delay. 


i 


& 2.  Preliminaries 


We  denote  w i t h *7 ( a , b ; R n ) , 1 < p <«>,  < a < b <«>,  the 

linear  space  of  functions  x:[a,b]  -»  Rn  such  that  |x|p  is 
integrable  on  [a,b].  |*|  denotes  the  Euclidean  norm  in 
Rn.  Lp(a,b;Rn)  is  the  Banach-space  of  equivalence  classes 
of  functions  in  <£p(a,b;Rn)  with  norm  II • II  . 

P D 1 0 C 

For  notational  purposes  we  also  introduce  the  spaces  LK’ 

Lp  ’ ^ oc  ( 0 ,oo;  Rn  ) of  equivalence  classes  of  functions  x:[0,«)  ->  Rn 
such  that  | x|  P is  locally  integrable,  1 < p<°°.  A family 

of  seminorms  on  Lp,^oc  is  given  by  ||x  ||  T = ( / 1 x I p dt)^p  , 

P 0 

T > 0.  C(a,b;Rn)  will  be  the  Banach-space  of  continuous  functions 

x:[a,b]  ->  Rn , -«>  < a < b < °°  , with  norm  ||x  ||r  = sup  I x I . 

L [a,b] 

It  will  not  be  necessary  to  indicate  the  interval  [a,b]  in  the 
various  norms. 


If  x is  a function  [-r,a) 


,a>0,0<r<°°,  then 


x^,  0 < t < a,  denotes  the  function  [-r,0] 


defined 


by  x.(s)  = x(t  + s),  s€[-r,0].X  will  be  the  linear  space 
Rnxi£p(-r,0;IRn),  1 < p < oo,and  X the  Banach-space  IRnxLp(-r,0;IRn) 
with  norm  ||(^,^)||  = ( |£>|  p+  II,/,  llp ) 1//p  . In  this  section  we  shall 
use  the  equ  i va  1 ent  norm  II  ) llQ  = 1^1  + lUIL  which  simplifies 
the  estimates  to  be  made  in  the  proofs.  'C  denotes  the  set 
{ (<p(0),<p)  l4>€C(-r,0;Rn)}  . 

Let  f be  a map^(f)  -*  Rn  where  i^c^f)  c and  take  g6Lp’^0C. 
Given  any  xQ  = (n,ip)€.2f  we  consider  the  Cauchy  problem 


x(t)  = f (x(t)  , x t ) + g(t) , t > 0 , 


(1.1) 


x(0)  = n,  x(s)  = ip(s)  a.e.  on  [-r,0). 


(1.2) 


A solution  of  this  problem  is  denoted  by  x(t)  = x ( t ; x"  , g ) or 


x(t)  = x(t;n><0»g)»  respectively,  and  is  a function 
[-r,°°)  -»  Rn,  a > 0,  such  that  (1.2)  holds,  x is  absolutely 
continuous  on  [0,a)  and  (1.1)  is  true  a.e.  on  [0,a). 

In  order  to  get  existence,  uniqueness  and  continuous  dependence 
of  solutions  we  impose  the  following  conditions  on  f: 

(HI)  If,  forsome  a>0,  x € L p ( - r ,a  ; Rn  ) and  x is  continuous 

on  [0,a]  then  the  map  t -♦  f(x(t),xt)  is  defined 

a.e.  on  [0,a]  , depends  on  the  equivalence  class 

1 n 

of  x only  and  is  in  L 1 ( 0 , a ; R n ) . 

(H2)  Given  3 > 0,  there  exists  a continuous  nonnegative 

and  nondecreasing  function  y(t)  = y^(t)  defined 
on  [0,<=°)  such  that  for  any  a > 0 the  inequality 

Jlf(x(s),x  ) - f(y(s),y  )lds  <y(t)(  / I x ( s ) -y ( s ) I pds  ) Vp 
o -r 

is  satisfied  for  all  t€[0,a]  and  all  functions 

x,y  in  Lp(-r,a;IRn)  which  are  absolutely  continuous 

on  [0,a]  with  ||x  II  < 3,  llyllp  <3- 

Remark  1 . 1 . Since  equation  (1.1)  should  incorporate  difference- 
differential  equations,  it  is  not  reasonable  to  demand  that 
c0(  f)  = 3f.  To  be  more  precise,  condition  (HI)  for  instance 
should  read  as  follows:  Given  x£L p ( - r ,a  ; R!n  ) such  that  there 
exists  a representation  of  x which  is  continuous 
on  [0,a]  then  for  each  representation  $ of  x the  map 

1 n 

t - f($(t),a  t ) is  in  the  same  equivalence  class  of  L 1 ( 0 , a ; Ht n ) . 

An  analogous  statement  is  appropriate  with  respect  to  condition 
(H2).  But,  since  expressions  like  f ( x ( s ) , x ) will  appear  only 
under  the  integral  sign,  the  formulation  of  conditions  (HI)  and 
(H2)  as  given  above  will  not  cause  any  confusion.  However, 
in  Section  2 we  shall  need  the  following  definition  of  f(n»<P) 
for  those  (n»ip)€X  such  that  there  is  a representation  (p  of  ip 
with  <p(0)  = n which  is  continuous  on  [ - r , 0 ] . In  this  case 


f ( n » 4> ) = f(ip(0),(p)  defines  f(n,<p)  uniquely,  because  tp  is 
uniquely  determined. 

Remark  1.2.  Let  equation  (1.1)  be  a difference-differential 
equation 

x(t)  = h ( x ( t ) , x(t-kj)  , . . . >x(t-km) ) , (1.3) 

ko=0  < k^  < ...  < km>  where  h is  defined  on  Rn(m+1)  and 
L i psch  i tz - co n t i nuou s , i.e.  for  each  M > 0 there  exists  a constant 
L(M)  > 0 such  that 

m 

1 h(yo’ ' • ' ,ym)_h(zo’  • * • ’zm)  1 5 L(M).Z  lyj~zj 1 

j =o  J J 

for  aii  (y0.---»ym)  and  (z0»--**zm)  such  that 


£ I y - 1 < M and  I |z.|  < M.  Suppose  there  exist  numbers 
j=o  J j=o  J 


p > 0 , a > 0 , 


and  q > 0 such  that 


L ( M ) < oM4  + p 


(1.4) 


for  M > 0.  Then 


J(t)  = J|f(x(s),xs)  - f (y(s) ,y  ) Ids 
m t° 

< I / (oiv  ( s ) +p  ) I x ( s - k . ) - y ( s - k . ) Ids, 

- j = o o J J 

mm  m 

where  6(s)=max(  I |x(s-k.),  L ly(s-k.)l)  < l|x(s-k.)|  + 

m d = 0 J J'  = 0 ..J 

+ l|y(s-k.)|.  Since  6£Lp(0,t;B!)  we  have  6 q ( s ) £L  1 " V p ( 0 , t ; R ) 

j =o  J 

provided  that  , 

q < p - 1 . 

Let  us  take  p = q+1.  Then  an  application  of  the  Holder  inequality 
gives 

|.  m t 

J(t)  < { P t1- Vp  + 0(J  6 p ( s ) d s ) 1 " 1/p  > I(/lx(s-k  .)-y(s-k.)  |pds)1/p 
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< (m+1  ){p  t1_  Vp  +o(/6p(s)ds)1'Vp}  ( Jlx(s)  - y(s)lpds)Vp  . 

o -r 

From  Minkowski's  inequality  we  get 

(/6p(s)ds)1/p  < I (/I  x(s-k  . ) I Pds).]/P  + 

O " j=0  0 J 

m t 1 < 

+ I (/ ly(s-k . ) I P) Vp  < 2 ( m+ 1 ) 3 

j = o o J 

i.e.  we  have  (H2)  where  y^  ( t ) is  given  by 

Yp(t)  = (m+1 ) (p  t1"  Vp  +o ( 2 ( m+1 ) 3 ) p” 1 ) . 

It  is  easy  to  see  that  also  condition  (HI)  holds. 


Remark  1.3.  Another  case  covered  by  conditions  (HI)  and  (H2) 
is  where  f is  a Li  pschi  tz-conti  nuous  map  Lp(-r,0;IRn)  -*  iRn,  i.e. 


If(<p)  " f(*)l  5 L (M)  Htp-«P  llp  (1.5) 

where  M = max  ( Ikpll  , Ityllp ) • Then  an  easy  calculation  shows 
that  (H2)  holds  with  y (t)  = tL(3). 

P 

Remar k 1.4 . Conditions  (HI)  and  (H2)  are  the  local  versions  of 
global  conditions  already  used  by  Borisovic  and  Turbabin  (6  1 
and  Banks  and  Burns  [2  ] in  the  linear  case.  The  conditions 
apply  to  functions  which  are  candidates  for  solutions  of  Cauchy 
problems  and  give  restrictions  for  the  right-hand  side  of 
equation  (1.1)  along  absolutely  continuous  prolongations  of 
initial  data.  For  existence  and  uniqueness  results  it  is  this 
behaviour  which  is  essential.  Related  existence  results  were 
obtained  in  [ 12]  and  [ 13] . 


Propos i ti on  1.1.  Suppose  that  (HI)  and  (H2)  hold  for  f and  take 
6 > 0.  Then  there  exists  a constant  a = a(6)  > 0 such  that 
the  Cauchy  problem  (1.1),  (1.2)  has  a unique  solution  x(t;x,g) 
existing  on  [-r,a]  provided  that  l|x,|l0  < 6 and  llg  llp  a 


< 6 . 
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Moreover,  we  have  I x ( t ) I < 6+1  on  [0,a]. 

Remark  1.5.  The  proof  of  this  proposition  will  show  that 
solutions  of  (1.1),  (1.2)  are  dependent  on  the  class  of 
x only.  Therefore,  it  is  legitimate  to  consider  initial  data 
in  X . 

Proof  of  Proposition  1.1.  Take  aQ  > 0 and  define  for  x = 
(n»<p)€X  with  llxll  < 6 the  function 


f o(s)  for  s€[-r,0) , 
f n for  s€ [0,aQ] . 


We  further  define 


S = {^6C(0,aQ  ;Rn)  I II# -n  llc  5 11  , 


where  n also  denotes  the  constant  function  on  [0,aQ], 

and  for  #€ S the  function 


?(s) 


f ip(s  ) for  se  [-r  ,0)  , 
| #(s)  for  s€[0,ao]  . 


The  operator  T is  defined  by 

t t 

(T*)(t)  = n+Jf (#(s)  ,#  )ds  + / g ( s ) ds , 
o o 


(1.6) 


t€[0,aQ],  #6S.  Note,  that  #^€1.  ^ ( - r , 0 ; Rn  ) for  s€[0,a  ]. 
Therefore,  according  to  (HI),  T#  exists  and  is  absolutely 
continuous  on  [0,aQ].  For  #eS  we  have  ||#||  < ( 6 + 1 ) ( 1 +aQ  ) ^ . 

Wechoose  y(t)  according  to  (H2)  for  3= ( 6+1 ) ( 1 +aQ ) . Then 
for  a^e(0,aQ]  and  #,x€S,  we  get  the  estimates 


k , 


l(T*)(t)-nl  5 y(«i)  (J1  l?(s)lpds)Vp  + 

- r 

+ /1|f(0,0)|ds  + / 1 1 g ( s ) I d s 
o o 

< Y(a1)3  + a1lf(0,0)|  + /i|g(s)|ds  = p^Oj) 

o 

and 

l(T*)(t)  - (TX  ) (t)  I < (J1 1 U.  ( s ) -x  ( s ) I pds ) 1/f5  Y(a,) 

o 

< «^p  y (a,  ) max  I *(s)~x(s)  I 
[O.aj] 

= p2(a1)  max  l*(s)-x(s)  I 
C [O.aj] 

for  alltetO.a,].  Since  Tim  p.(a)  = 0,  i = 1,2,  these  estimates 
1 a-»0+  1 

show  that,  for  sufficiently  small,  (1.6)  defines  a contraction 
S0  ->■  So,  So  being  the  complete  metric  space 

So  = {^ecco.ap*")  I llip-n  llc  < 1}  . 

Therefore,!  has  a unique  fixed  point  in  SQ  which  gives  the 
unique  solution  to  the  Cauchy  problem  (1.1),  (1.2)  with 
llxll  <&  Existence  of  these  solutions  is  guaranteed  on  each 
interval  [0,a^]  where  is  determined  by  the  conditions 
p 1 ( a 1 ) - * and  pp(ai)  < definition  of  SQ  it  is  clear 

that  |x(t;x,g)|  < 6 + 1 on  [ 0 , a j ] . 

Remark  I.6..  The  proof  of  Proposition  1.1  shows  that  in  order  to 
get  existence  and  uniqueness  of  solutions  to  the  Cauchy  problem 
(1.1),  (1.2)  condition  (H2)  can  be  weakened.  The  estimate  in 
(H2)  only  needs  to  hold  for  those  functions  x,y  which  satisfy 
in  addition  x(t)=y(t),  t€[-r,0).  For  instance  this  weaker  version 
of  (H2)  holds  for  the  scalar  equation 

x(t)  = sgn  x(t-l ) 


whereas  (H2)  does  not  hold.  The  full  power  of  condition  (H2) 
is  needed  to  prove  continuous  dependence  of  solutions  on 
initial  data.  It  is  easily  seen  that  in  case  of  the  example 
given  above  solutions  do  not  depend  continuously  on  inital  data. 

Remark  1.7.  If  condition  (H2)  is  a global  condition,  i.e.  the 
function  y(t)  does  not  depend  on  3,  then  solutions  of  (1.1),  (1.2) 
exist  globally.  This  is  true  because  we  can  take  S = C(0,ao;Kn) 
and  p^fa)  is  only  dependent  on  a and  y(t).  This  implies 
that  each  solution  can  be  prolongated  on  an  interval  of  length 
o<j  where  is  determined  by  P2(a^)  < 1. 

Propos i ti on  1.2.  Suppose  that  cond i ti ons  ( HI ) and  (H2)  hold. 

a)  If  x(t)  = x(t;n»4>»g)  and  y(t)  = x(t;n,<P,g)  are  two  solutions 
existing  on  the  same  interval  [-r,T]  then  there  exist 
nonnegative  constants  k Q , k ^ and  a such  that 

I x ( t ) -y  ( t ) | < K0ll(ri,(p)-(~,^)ll0eot+K1llg-g||p)Teat  (1.7) 

for  t£  [ 0 , T ] . The  constants  are  dependent  on  T , llx  ( • ) ||  y 
and  lly(-)Hp- 

b)  Let  x ( t ;n  »ip»  g ) be  a solution  with  maximal  interval  of  existence 
[-r,t+) , 0<t+  < °°.  Then  for  every  Te(0,t+)  there  exists  a 
positive  constant  xQ  such  that  x(t;~,(p,g')  exists  on 

[-r,T]  provided  that  ||(~  ,£)  - ( n ,<P)  IIQ+  llg-g  llp  ^ T < tq.  The 

constant  t is  dependent  on  sup  ||(  x ( t ) , xt ) II  . 

0 [0 , T ] 1 0 

c)  Let  x(t;n»<P.g)  be  as  in  b)  and  suppose  that  (nn.<Pn)  and  (g  ) 

i 1 1 II  II 

are  sequences  in  X and  L*5’  oc  such  that  (nn>tpn)  -»  (n.»p)  and 

||g  -g  ||  T -*  0 for  each  T > 0.  Then 
n p 9 I 

lim  x(t  ;nn,<Pn»gn)  = x(t;n,«P,g) 

n-K» 

uniformly  on  compact  subi ntervals of  [0,t+). 


Proof . We  choose  3 > 0 such  that  ||x  11^  < 3 and  ||y  ||p  <3- 
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Then  with  y(t)  = Y^(t)  according  to  (H2)  we  obtain 


I x ( t ) - y (t ) I < ln-n!  + y (T)  Ikp-Jpll  + y (T ) (/ I x ( s ) -y  ( s ) I pds  ) ^ 

p o 


+ T1"^ 


iig-giir 


If  we  define  v(t)  = |x(t)-y(t)|, 
t 

V ( t ) = Jvp(s)ds  for  t > 0 and  p = In-rTl  + y(T)  lltp-cp II 
o " • 

+ T 1 " ^P  Ilg-gIL  t we  get 

v(t)  < P + y (T)  (V(t))1^  , t£[0,T]  , 


V(t)  = vP(t)  < (p+y(T)(V(t))Vp)P,  te [0,T]  . 


(1.8) 


In  case  p=l  Gronwall's  inequality  gives  the  desired  result. 

In  the  general  case  we  have  to  proceed  as  follows.  Substitute 
y(t)=y(T)V(t)1/p  . Then  integrating  (1.8)  we  obtain 


y(T) 


’ V 

( (y(Myjp  1[ly  - 1 ■ te[0-Tl- 


If  we  note  that 

(p+y)p 


00  „ i / i V 

! (P-M"1)  P 

~ v / , ,v+l 
v=o  (y+p ) 


for  nonnegative  y then  termwise  integration  gives 


ln(P+Y(T)VVp  )+  ” 


< (y(T)/P)t  + In  p+  z 

V = 1 

From  this  we  get  the  estimate 


OO  / 1 \ V T 1 

v(-l)  (P-1W  P 

7 V { V ' {~  7T7p  } 

v = l p +y  j V r 

00  / i \V  + 1 „ , 

In  ri  4.  T ■(  ~ • LI- /P”^\ 


r )• 

' V ' 


ln(P+y(T)V1/p  ) < (y(T)p/p)t  * 1 np  +2  7 ^(P'^) 

k = o 


v ( t ) < p +y  ( T ) V 


< p~  exp(y(T)P/p)t  ,te[0,T], 
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l where  ~ = exp(2  I 2~k+I  ( 2k  + l ^ • Note,  that  the  series 

k = o^ 

appearing  in  this  proof  are  convergent.  According  to  the  defini- 
tion of  p and  v(t)  the  last  estimate  is  just  (1.7)  with 
o=y(T)^/p,  kq  = ic  max  ( 1 , y (T ) ) and  = < T . 

In  order  to  prove  part  b)  we  choose  6 > 0 such  that 

||(  x ( t ) , x . ) II  < 6 on  [0,T]  and  ||g  II  T < 6 . Then  according  to 

Proposition  1.1  there  exists  a constant  a > 0 dependent 

on  6 only  such  that  each  solution  y(t)  = x ( t ; ~ , (p , g” ) exists  on 

[-r,a]  provided  that  ll(~,(p)  II  < 6 and  llg  II  „ < llg  lln  T < 6 . 

Moreover  I y ( t ) I < 6 + 1 on  [0,a].  Then  by  part  a)  of  Proposition 
1.2  there  exists  a constant  t ^ > 0 (dependent  on  6 only) 
such  that 

ll(n,tp)  - (~,£)  llQ  + Hg-gll  T < tj  (1.9) 

implies  ll(y  (a)  ,y  ) llQ  < 6.  Then  by  Proposition  1.1  all  solutions 
y(t)  with  (1.9)  exist  on  [-r,2a]  with  ly(t)|  <6  + 1 on  [0,2a]. 

By  part  a)  of  Proposition  1.2  we  can  choose  a < tj  such 
that  Il(y(2a)  ,y2a)  llQ  < for  all  solutions  y(t)  such  that  (1.9) 
holds  with  replaced  by  t Proceeding  in  this  way  we  arrive 
in  a finite  number  of  steps  at  a number  tq  > 0 which  has  the 
properties  stated  in  part  b)  of  the  proposition. 

f 

Part  c)  is  an  easy  consequence  of  parts  a)  and  b). 

Remark  1.8 . If  condition  (H2)  is  a global  condition  then  an 
inspection  of  the  proof  given  for  part  a)  shows  that  the  constants 
a,  kq  and  are  independent  of  x(*)  and  y(*). 
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I & 3.  The  nonlinear  local  semigroup 

In  this  section  we  consider  the  autonomous  equation 

x(t)  = f(x(t),xt),  t > 0.  (2.1) 

Conditions  (HI)  and  (H2)  are  assumed  to  hold  for  f throughout 
this  section.  Civen  xo  = (n,<p)€X  we  simply  write  x(t;n»<p) 
or  x(t;,x’  ) instead  of  x(t;n»<P,0)  or  x(t;xo,0),  respectively. 

To  each  x€X  there  corresponds  the  maximal  interval  of  existence 
t-r,t~)  for  x(t;x)  where  0<t~  < °°.  It  is  clear  that 
(x(t;x) ,xt(x) )€X  for  all  t€[0,t~). 

Therefore  (t,x)  -*  (x(t;x)  ,Xj.(x;) ) defines  a map 
T:D  -*  X where  D = { ( t , x ) I x€X , t€ [ 0 , t~) } . 

Propos i ti on  2.1.  v i s a local  semidynamical  system,  i.e. 

: (i)  f(0,x)  = x for  all  x£X, 

(ii)  4'(t^  + t2;x)  = y ( t^  ;y  ( tg  ;x") ) for  all  x€X  and  all  nonnegative 

t ^ , 1 2 such  that  tj+t2  < t~, 

(iii)  v is  continuous  on  D, 

' (iv)  for  any  sequence  (>T  ) such  that  7 -»  x€X  we  have 

liminf  t~  > t~ 
n-*~  n 

Proof . (i)  is  evident  by  definition  of  f and  (ii)  follows  imme- 
diately, since  equation  (2.1)  is  autonomous,  (iv)  is  a conse- 
quence of  Proposition  1.2,  part  b).  In  order  to  prove  (iii) 
let  ( t , x ) , n = 1 , 2 , . . . , be  a sequence  in  D with 
( t n , x”n ) -♦  (t,x)eD.  Since  tn  -»  te[0,t~)  it  is  clear  that  there 
exists  a number  nQ  such  that  (tn,x)6D  for  all  n > n . 

Proposition  1.2,  part  c)  implies  IlY  (tp  ;xp  )-¥(  tn  ;*?)  ||  0 as 

n -♦  oo.  Therefore  it  is  sufficient  to  show  that 

lim  y ( t n ; ST)  = H’(t;x).  But  this  is  clear  by  continuity  of 
n-*oo 
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x(t;x)  for  t > 0 and  by  II x t (x)-xt(x)||  -*  0 as  n 

We  shall  deal  with  the  operStors 

Tt  = v(t;.  ) , t£[0,t  + ) , 

where  t+  = sup  t-i.  If  D.  denotes  the  domain  of  T*  then  it  is 

xex  / 1 1 

clear  that 

Dt  = (x€X| (t,x)6D)  , t€  [ 0 , t+ ) . 

We  list  the  following  properties  of  the  family  { T 1 1 1€ [ 0 , t+ ) } : 

c for  0 < t^  < t^  < t+,  Dq  = X (2.2) 

and  x€Dt  for  all  te[0,t~). 

For  all  t > 0 the  domain  Dt  is  an  open  subset  of  (2.3) 

X and  for  each  6 > 0 there  exists  a t >0  such  that 

o 

{XI  lixll  < 6}  c Dt  for  all  t€[0,t  ]. 


T.  ..  x = T.  ( T . x)  for  all 

M + t2  C1  +r2 

that  t ^ + 1 2 < t and  all  x£D 
For  t€[0,t+)  the  operator  T 


nonnegative  tj,t£  such 


t is  continuous  on  D^. 


The  map  t -*  Tt  is  strongly  continuous  on  [0,t+), 

i.e.  t -»  T.x  is  continuous  for  all  x€X  on  [0,t~)and 

^ y 

Tq  = I (=  identity  map  on  X) 


(2.4) 


(2.5) 

(2.6) 


Properties  (2.2),  (2.4),  (2.5)  and  (2.6)  are  trivial 
consequences  of  properties  proved  for  v . Property  (2.3)  is 
clear  from  Proposition  1.2,  part  b)  and  from  Proposition  1.1. 


Remark  2.1.  Regarding  continuity  of  Tt  on  we  already  have 
proved  more  than  just  continuity.  Combining  Proposition  1.1 
and  Proposition  1.2,  part  a),  we  see  that  for  each  6 there 
exists  a t > 0 such  that  a Lipschitz  condition  holds  for 
on  { ST | ||x|lo  < 6}  provided  that  te[0,t  ].  Moreover,  a 
combination  of  parts  a)  and  b)  of  Proposition  1.2  shows  that 


L 
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for  each  xGX  and  each  T€[0,t~)  there  exists  a neighborhood 
U of  x such  that  a Lipschitz  condition  holds  for  on  U 
provided  that  t € [ 0 , T ] . From  this  by  standard  compactness 
arguments  we  obtain  that  Tt  is  Lipschitz  on  compact  subsets 

Remark  2.2.  If  (H2)  is  a global  condition  then  the  operators 
Tt  are  defined  on  all  of  X for  all  t > 0 (c.f.  Remark  1.7) 
and  {Ttlt  >0}  is  a strongly  continuous  semigroup  of  continuous 
operators  . 

We  call  (Tt I te[0,t+)}  the  local  semigroup  associated  with 
equation  (2.1).  The  infinitesimal  generator  A of  this  semi- 
group is  defined  by 

Ax  = lim  |(T.x  - x)  , xe^(A),  (2.7) 

t-»0+  r z 

where  c#(A)  is  the  set  of  all  x'eX  such  that  the  limit 
in  (2.7)  exists. 

We  shall  need  the  following  additional  condition  on  f: 

(H3)  The  map  t->  f(x(x),x  ) is  continuous  on  [0,a)  for  all 
functions  x(-)  which  are  absolutely  continuous  on 
[-r,a),  where  a is  any  positive  number. 

Proposition  2.2.  Suppose  that  (HI)  - (H3)  hold. 

Then 

a)  x = (n»<p)€X  is  in^(A)  if  and  only  if 

(i)  <peW1,p  (-r ,0;Rn) , i.e.  <p  is  absolutely  continuous 
on  t-r,0]and  tpG Lp,(-r,0;R5n)  , 

( i i ) ip(  0 ) = n . 

b)  Ax  = (f  (tp(0)  ,(p)  ,(p)  for  all  x = («p(0)  ,<p)e«0(A)  . 

Proof . Suppose  that  x = (n.<p)  is  in  t3(A)  and  put  A(n,<p)  = 

(^ ,<i> ).  This  is  equivalent  to 


lim  i(x(t;x)-n)  = > in  Rn 
t-0+  1 


(2.8) 


(2.9) 


lim  T(xt(x)-4>)  = 'P  in  LH(-r,0;R  ). 
t-0+  1 1 

Relation  (2.9)  is  equivalent  to 
o , 

lim  / I f-(xt  (x  ) ( s )-<p(  s ) )-i|/  ( s ) I pds  = 0 
t-»0  + -r  1 1 

or  to 

-t  , 

lim  / I 4-(tp(  t + s ) -4>(  s ) ) -ij;  ( s ) I pds  = 0 (2.10) 

t-*0  + -r  1 

and 

o , t+s 

lim  / 1 4- ( n + f f ( x (t  ; x ) , x (x))dT  - cp(s))-ip(s)|Pds  = 0.  (2.11) 

t-0  + -t  1 o T 

i 

By  a characterization  of  functions  in  W 1 ’ p ( - r , 0 ; Rn ) 

(of.  [5  ] , p.21,  or  [7  ]»  p.  154)  relation  (2.10)  is 
equivalent  to  statement  (i)  in  the  proposition  with  <p  = ^ . 

Relation  (2.11)  means  that  for  each  e > 0 there  exists  a 
t > 0 such  that 

0 t+s  ,, 

( / I n + / f (x  (t  ;x)  ,x  (x)  )d  x-tp(  s )-t^(  s ) I pds  ) vp  < ct  (2.12) 

-t  o T 

for  all  te [0, t ] . 

Suppose  that  4)  is  continuous  and  put  6 = I cp  ( 0 ) - n I . Then 

for  t sufficiently  small  we  have  the  estimates 

o t+s  ,, 

( Jln+  / f (x(t  ;x)  ,x  (x)  )di-4>(s)-t^(s)  I pds)i/p  (2.13) 

-to  T Q 

< 26  t1/p  + t ( / U ( s ) I pds  ) 1/p 
-t 

and 

o t + s _ 1/n 

( / I n+  / f(x(t;x),x  (5f))dr  - 4>(s)  - t<|»  ( s ) I pds  ) Vp  (2.14) 

-t  0 T 

> -ifit1^  - t ( /U(s)  |pds) Vp  . 

" -t 

If  6 = 0 then  (2.13)  shows  that  (2.12)  is  true.  But  if 


^ — 1 


<5  > 0 then  by  (2.14)  we  see  that  (2.12)  cannot  hold. 

So  (2.13)  is  equivalent  to  ip(0)  = n provided  that  ip  is 
continuous. 

Since  x(t;x)  = n + /f(x(t;x),x  (x))dT  for  t > 0,  relation 

o T 

(2.8)  is  equivalent  to 

& = f (<p(  0)  ,cp) 

provided  that  ip  is  absolutely  continuous  and  tp ( 0 ) = n . Thus, 
the  proof  of  Proposition  2.2  is  finished. 

Remark  2.2.  Each  paper  dealing  with  semigroup  theory  and  functio- 
nal-differential equations  in  the  state  space  X contains  a 
result  or  a definition  corresponding  to  Proposition  2.2  (cf. 
for  instance  [16],  [20]).  The  proof  is  quite  standard 

and  only  included  in  order  to  make  this  paper  more  sel f conta i ned . 

Remark  2.3.  The  additional  condition  (H3)  is  not  very  stringend. 
In  fact  all  important  types  of  retarded  functional -di fferential 
equations  are  incorporated.  If  we  only  are  interested  to  prove 
the  result  in  Proposition  2.2  it  is  sufficient  that 
x -»  f (x(t  ) ,x  ) is  right-hand  continuous  at  t = 0 for  the  class 
of  functions  appearing  in  this  additional  condition.  But  the 
stronger  form  of  this  additional  condition  is  needed  in  order 
to  prove  the  next  result. 

Proposition  2.3.  Suppose  that  (HI)  - (H3)  hold.  Then 
a)  TtxG(D(A)  for  all  x€^(A)  and  all  t€[0,t~). 
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Te(0,t~) . The  proof  of  part  a)  is  finished  if  we  note  that 
( ip(t  + s)  for  t + s < 0, 

*t(x)(s)  = j 

f ( x ( t + s ;x  ) , xt  + s ( x ) ) for  t + s > 0. 

In  order  to  prove  part  b)  we  write 

ATtx  = (f(x(t;x),xt(x)),xt(x)) 

By  (H3)  it  is  clear  that  the  first  component  is  continuous 
with  respect  to  t.  Continuity  of  x^x”)  as  a map 
[0,t~)  -*  L p ( - r , 0 ; R!n  ) follows  from  x ( • )€L  p ( - r , T •,  IRn  ) for  all 

T6[0,t~)  and  the  continuity  the  LH-norm  with  respect  to  trans- 
lation. 

In  order  to  proove  part  c)  we  first  have  to  show  that  the  map 
t -*  Tj.x’  is  differentiable  a.e.  on  (0,t~)  provided  that  x£a5(A). 

This  is  trivial  for  the  first'  component  x(t;x”)  of  Ttx\  With  respect 
to  the  second  component  Xj.(x')  one  has  to  show  that 

lfH(xt  + h(")'xt(3r^  ‘ xt(x)"p  = °- 

We  leave  this  to  the  reader.  Again,  one  has  to  use  the  result 
given  in  [7  ] , p.  154. 

Remark  2.4.  Proposition  2.3  shows  that  with  equation  (2.1)  we 
can  associate  an  abstract  evolution  equation  in  X.  Moreover, 

Tt?  is  a strong  solution  of  this  evolution  equation  for  x£o)(A). 

2)(A)  is  a linear  subspace  of  X not  dependent  on  f.  It  can  be 
characterized  as  that  set  of  initial  data  x corresponding 
to  solutions  x(t;x”)  which  are  absolutely  continuous  with 
Lp- i ntegrabl e derivative  on  compact  subintervals  of  t-r,t~). 

If  we  take  the  state  space  C(-r,0;Rn)  instead  of  X then  <3(A) 
is  dependent  on  f and  is  given  by  {ipeC(-r,0;Rn)|ip£C(-r,0;ll!n) 
and  <i>(0)  = f(ip(0),ip)}  whereas  Aip  = tp.  Here,  S(A)  can  be 


characterized  as  the  set  of  those  initial  data  x corresponding 
to  solutions  x(t;x’)  which  have  a continuous  derivative  on 
[-r,t~)  (cf.  for  instance  [19]). 

Proposition  2.2  shows  that  the  operator  A can  be  written  as 


where 


A = A + A, 
o i 


Ao(<P(0),tp)  = (0,ip) 


(2.15) 


AjCipfOJ.tp)  = (f  (<p(0)  ,<p)  .0) 

for  all  (tp(  0 ) ,ip)€0(  A ) . Of  course,  A1  is  defined  on  all  of 
£=  {(<P(0),<p)|<p€C(-r,0;IRn)^. 

'C  will  be  considered  as  a Banach-space  with  norm  II ( cp ( 0 ) ,cp)  ||  = 

lkPllc. 

Aq  is  the  infinitesimal  generator  of  the  linear  semigroup 
[Stlt  > 0$  generated  by  the  solutions  of  x = 0,  i.e.  St(n,ip)  = 

( n > ^ ) where 

r n for  t+s  > 0, 

*(s)  =< 

( ip(  t + s ) for  t + s < 0 . 

Since  X - £){&)  =«?(AQ)  we  have  from  linar  semigroup  theory 
(see  for  instance  [21]) 

Proposition  2 . 4 . is  dense  in  X. 

Splitting  up  the  infinitesimal  generator  A as  in  (2.15)  shows 
that  properties  of  f in  a straightforward  way  lead  to  properties 
of  A^  and  therefore  of  A.  This  is  one  advantage  of  taking  X as 
the  state  space.  As  we  have  stated  above^(A)  is  dependent  on  the 
right-hand  side  of  equation  (2.1)  and  properties  of  f are  re- 
flected intv(A)  if  we  take  the  state  space  C ( - r , 0 ; !Hn  ) . 
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& 4.  Representation  of  the  solution  semigroup,  global  case 

Throughout  this  section  we  assume  that  (HI)  - (H3)  hold  for 
f and  that  (H2)  is  a global  condition.  We  shall  make  use  of 
some  results  of  nonlinear  semigroup  theory  (as  a reference 
see  for  instance  [5  ] ) . 

* 

If  Y is  a Banach  space  and  Y is  its  dual  then  the  duality  map 
F on  Y is  defined  by 

F ( x ) = {x*€Y*|  llxll2  = llx*  II2  = <x,x*>] 

for  x£Y,  where  denotes  the  pairing  between  elements  of 

Y and  Y . If  Y is  strictly  convex  then  F is  univalent.  An  ope- 
rator AcYxY  is  called  dissipative  if  <u1-u2,x*)>  < 0 for  all  pairs 
[Xj.u^GA,  i = l ,2,  and  some  x*£F(x1-x2)  . If  A and  F are  univalent, 
the  dissipativeness  condition  is 

“CAx^  - Ax2  , F(Xj-x2)>  < 0 (3.1) 

for  all  x i € ©5(A)  , i=l,2. 

An  equivalent  cha rac ter i za t i o n of  an  dissipative  operator  was 
given  by  Kato  [14]:  A-cal  is  dissipative  if  and  only  if 

IIXj-AUj  -(x2-Xu2)ll  > ( 1 -coX  ) llx  j - x g II  (3-2) 

for  all  [ x .j , u .j  ]€  A and  all  A€(0,~).  If  (I-AA)'1  exists  on 
Y,  then  (3.2)  implies 

|[(I-AA)"1x1  - ( I -AA  )~  1x2  II  < -j^—IIXj-x^  11 

for  x.£Y,  i = 1 , 2 , and  \€(0,-^). 

The  spaces  X supplied  with  the  norm  ll(i°!,4’)ll  = ( l/th  IliHIp)  ^ 

and  X*  = RnxLq(-r,0;Rn) , = 1,  supplied  with  the  dual  norm 

are  uniformly  convex.  Therefore  the  duality  map  on  X is  univalent. 
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In  [ 20] , Proposition  2.1,  it  was  proved  that  the  duality  map 
F ( x ) for  x = (n,<p)€X  is  given  by 

<(*>.*). F(x)>  (3.3) 

= llx||2“P(^TnlnlP_2+  J ip(s)T4>(s)  lcp(s)  1 p_2ds) 

- r 

for  all  (A^  )£X. 

With  respect  to  the  linear  operator  Aq  we  have  the  following 
results. 

Propos i t i on  3 . 1 . 

a)  ( I -x AQ ) ” exists  on  X for  all  x > 0 and  is  given  by 

(i-xa  rV,*)  = (Ae*/x^4/e‘T/x^T)dT))  (3-4) 

0 Ao 

for  all  (s?,ip)€X  and  x > 0. 

b)  AQ-|l  is  dissipative  and  therefore 

||(I-xA0)‘1x  ||  < Mx  II 

obtains  for  all  xeX  and  x£(0,p). 

c)  The  restriction  of  (I-XA0)_1  to  the  Banach  space  f i s a 
contraction,  i . e . 

||(I-xAo)'1x|lif  < llxll^ 

for  all  'x  € £. 


Proof . The  semigroup  S , t > 0,  with  infinitesimal  generator 
Aq  is  equ icont  i nuous , i.e.  for  some  constant  M > 0 we  have 
||S  t II  < M for  all  t > 0.  Therefore  existence  of  (I-xAQ)  * for 
all  x > 0 immediately  follows  from  linear  semigroup  theory. 
The  representation  (3.4)  is  obtained  by  solving  the  equation 
( I -x  A ) (n.ip)  = (Ai|>)»  (A,  ip ) GX , which  is  equivalent  to  n = 
and  ip(s)  - x<p(s)  ip(s)  , s€[-r,0]. 


With  respect  to  part  b)  we  refer  to  Proposition  3.1  of  [20]. 
Part  c)  immediately  follows  from  the  fact  that  the  restriction 
of  the  semigroup  , t > 0,  to  the  space  'C  i s a contraction 
semigroup.  But  the  result  can  also  easily  be  obtained  by  direct 
computation  using  (3.4). 

Remark  3.1.  From  Proposition  3.1,  part  a),  we  infer  that  Aq  is 
m-dissipative,  i.e  ^(I-aAo)  = X for  all  A > 0. 

As  stated  above  we  deal  with  the  case  where  (H2)  is  a global 
condition,  i.e.  the  function  y ( t ) is  not  dependent  on  3-  For 
e > 0 we  define  the  operator  A with  ^(A  ) = 3)  by 

A.(«p(0),«p)  = (7/f  (S  (<p(0)  ,(p)  )ds  , <i>)  (3.5 

e e0 

A corresponds  to  equation 
1 e 

x(t)  = 7 Jf  ( S s ( x ( t ) ,xt))ds.  = fe(x(t),xt)  (3.6 

in  the  same  way  as  A corresponds  to  equation  (1.1). 

Note,  that  f£  is  defined  on  all  of  X,  if  for  (n,<P)€x  we  put 

f (n>4>)  = — Jf  ( S (n  ,«p)  )ds . If  (H3)  holds  for  f then  it  is 

also  true  for  f . 

c 

Proposition  3.2.  Suppose,  that  the  function  y(t)  in  (H2) 
does  not  depend  on  3. 

a)  The  operator  A e I is  dissipative  witho=  ( e ) » 

0 < e < 1,  and  cj  = ^ *y(0  . e>l. 

b)  The  solutions  x^e^(t;  x)  of  equation  (3.6)  exist  globally 
for  all  xeX  and  ( x ^ e ^ ( t ;x ) , x ^ ( x”) ) = T^x  where 

T^x  = 1 im  U-7Ae)~nx  . (3.7 

n-*«> 

the  limit  existing  uniformly  on  compact  intervals. 

c)  t|£^x  is  the  unique  solution  of 
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T<c>x  = Stx  + {S 


t-s\,lTsE>*  ds  • 1 i 0 • 


where  'xeX  and  A , = A 

C j i £ 


V 


(3.8) 


Proof.  Choose  7.  = (<Pi  (0)  ,<p . )€<0  , i -1 , 2 . Then  by  (3.3)  we  get 

<Al,t  *1  ' Al,t*2  ' F(xrx2)> 

= |llx1-x2ll2"pHP1(0)-<(.2(0)lp'2(®)(0)-t(.2(0))T/(f(Ss(<p1(0),«.1) 


- T(Ss(cp2(0)  ,cp2)  )ds 

< | Y (e  ) ll'x1-x2  ll(elu>1(0)-cp2(0)|p+  / I <Pj  ( s ) -<p2  ( s ) I Pd  s ) 1/p 


|Y(e)||xrx2ll2  for  0 < e < 1, 


-1y  (c  ) llxrx2 


for  e > 1. 


This  estimate  and  the  dissipativeness  of  AQ-il  implies  the 
desired  result. 

Since  A£-u>I  isdissipative,  an  application  of  Theorem  I in 
[9  ] shows  that  A£  is  the  infinitesimal  generator  of  the 
semigroup  given  by  (3.7)  provided  that<5£((I-  A£)  ^)  = 


X 


for  all  sufficiently  small  x.  In  order  to  verify 


this  we  have  to  discuss  the  equation  x = (I-xA£)y  , yeJ)  , xeX, 
which  is  equivalent  to  the  fixed  point  equation 

y = (I-XAJ1?  + A(I-xAo)_1Aey 

for  y.  It  is  easy  to  see  that  this  fixed  point  equation  has  a 
unique  solution  ye^for  all  x€X  provided  thatxe(0,^).  It  was 


shown  in  [20),  Proposition  5.9,  that  (x^e^(t;x)  , x£0/(x))  = 
T^x  for  all  xeX  and  t > 0. 

Formula  (3.7)  also  follows  from  the  results  given  by  Webb  in 
[18).  There  it  is  also  proved  that  t[e^x  is  the  unique  solution 
of  (3.8)  . 


(O 
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Remark  3.2.  It  is  immediate  that 

lim  Ax:  = Ax;  for  all  ”x  e . 
e-»0+ 

If  lim  -y(e)  = y'(0)  exists  then  from  the  proof  of  Proposition 
c-O-i  e 

3.2  we  obtain 

<Axx-Ax2  , F(x1-x2))>  < (|+  y ' (0) ) llxj-x^l2 

for  all  *.  $-T)  , i = 1 , 2 , i.e.  A-talVii  ssi  pati  ve  with  u=i+y'(0). 

But  in  general  y'(0)  does  not  exist.  The  only  case  where  y'(0) 
exists  seems  to  obtain  when  condition  (1.5)  holds  for  f 
(with  L independent  of  M ) . In  this  case  y ' ( 0 ) = L.  But  then  an 
easy  calculation  directly  shows  that  A-(-i+L)I  is  dissipative. 

In  case  of  the  difference-differential  equation  (1.3)  y'(0)  does 
not  exist,  for  instance. 

Remark  3.3 . If  the  right-hand  side  of  equation  (2.1)  is  of 
special  form,  then  it  was  shown  in  [20]  that  A -cal  is  dissipative 
for  some  u > 0 provided  one  uses  weighted  norms.  The  weighting 
function  then  is  defined  by  the  right-hand  side  of  equation 
(2.1). 

Remark  3.4.  In  order  to  prove  existence  of  solutions  for  equation 

(3.6)  it  is  not  necessary  to  use  the  Crandal 1 -Liggett-Theorem 

and  Proposition  5.9  of  [20].  From  the  properties  of  f it  follows 

that  f is  defined  on  all  of  X and  is  globally  Lipschitz  with 
e 1 

Lipschitz  constant  — y(e)  for  0 < e < 1.  By  direct  verification 
one  can  show  that  also  conditions  (HI)  and  (H2)  hold  for 
equation  (3.6).  A corresponding  function  y(t)  in  (H2)  is 
Y(e)(t)  = iY(e)t(l+(f)Vp  ),  0 < £ < 1. 

Thus,  existence  and  uniqueness  of  solutions  to  equation  (3.6) 
follow  from  the  results  of  Section  1. 

But  the  functions  y ^ e ^ ( t ) , are  dependent  on  e and  cannot  be 
estimated  by  an  function  independent  on  e as  e -*  0 . 
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ir  ) 

In  order  to  prove  that  the  semigroups  T£  ' approximate  the 
solution  semigroup  Tt  we  need  additional  conditions: 

(H4)  f is  locally  Lipschitz  on  £f,  i.e.  for  all  3>0 
there  exists  a constant  L(3)  > 0 such  that 

I f (tp  ( 0 ) ,cp ) - (0)  ) | < L ( 3 ) - lltp-ip  Un- 

provided that  liipll^  < 3 and  ||i{;  IIq  5 3- 

.( H 5 ) For  all  3 > 0 there  exists  a function  Yp(t)  which  can 
be  used  in  condition  (H2)  not  only  for  equation  (2.1) 
but  also  for  equation  (3.6)  if  e is  sufficiently  small. 

Remark  3.5 . Condition  (H4)  obviously  implies  (H3).  In  case  of 
equation  (1.3)  conditions  (H4)  and  (H5)  are  true  without  any 
further  assumption  on  h.  This  is  clear  with  respect  to  (H4). 

In  order  to  prove  that  (H5)  is  true  one  has  first  to  inter- 
change the  order  of  integration  and  then  to  proceed  in  the 
same  way  as  in  Remark  1.2.  We  see  that  the  function  y^(t)  given 
there  can  be  used  in  (H5)  for  e€(0,k^]. 

If  (1.5)  holds  for  f then  again  (H4)  and  (H5)  are  true  without 
further  restrictions  for  f. 

Remark  3.6.  If  condition  (H4)  is  true  for  f then  it  is  also 
true  for  the  equations  (3.6)  with  the  same  Lipschitz  constant 
L(3).  This  follows  from  the  estimate 

lfe(4>(0),<p)-fe(*(0),*)|  < | / I f ( S s ( <p  ( 0 ) »cp)-f  (S$(ip(0)  » tp ) ) Ids 

< 7L  ( 3 ) / IIS  (ip(  0 ) - ij;  ( 0 ) , tp- ip ) ||  ds 

e 0 t 

< L ( 3)  lltp- llC)B  = max(||ip|lc,  ||<|;  11^),  where  we  have  used  the 
fact  that  St  restricted  to  )f  is  a contraction. 


Propos i tion  3.3 . Suppose  that  conditions  (HI)  - (H5)  hold  for  f 
and  that  the  function  y(t)  appearing  in  (H5)  can  be  chosen 
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independent  of  3.  Then 

1 im  T |E  ^ x = T . x 
e-*0+  E E 

for  all  xGX,  the  limit  being  uniform  on  compact  intervals 


Proof . It  is  sufficient  to  prove  lim  x^  '(t)=x(t)  uniformly 

e-*0+  /\  / \ 

on  compact  intervals,  where  x(t)=x(t;<x)  and  x ' E ' ( t ) =x ' E ’ ( t ; x ) 
are  the  solutions  of  equation  (2.1)  and  (3.6)  with  initial 
data  3c,  respecti  vely. 

We  first  consider  the  case  3c£  'C  . By  assumption  (H4)  we  have 

|/f(Ss(x(e)(T),XT(E)))ds  = f(Ss+(x(c)(x),xT(e)))  (3.9) 

afc  jfc 

for  t > 0 where  s = s (x)G(0,e). 

We  see  from  Remark  3.6  that  the  set  { x ^ e ^ ( - ) I e > 0}  is  an 
equ i co n t i nuous  set  of  functionson  each  interval  t-r,T], 

T > 0.  Therefore 

J||(x(e)(x),x.|e))  - Ss+(x(E)(x),x|e))||  dx  < h(c;t)  (3.10) 

0 ^ 

for  e > 0,  t > 0 , wiiere  lim  h(e; t)  = 0. 

e-»0  + 

From  (3.9)  and  (3.10)  we  get  the  estimate 

1 X ( t ) - x(e)(t)l  < / I f (X(x  ) ,x  ) - f(s s + (x(e)(x),x|e)))ldT 

0 

5 Jlf(x(T),x  ) - f(x(E)(x),xJE))|dx 
0 

+ /|f(x(E)(x),x(E))  - f(Ss*(x(E)(x),x.[E)))|di 
0 

t , . 

< L(«)J  II  x - x<e  > llpdx  + L ( 6 ) h ( e ;T  ) 
o T 1 L 

for  te[0,T),T  > 0,  where  6 is  a uni  form  bound  for  x(-)  and 
x^E  ^ ( • ) , e > 0,  on  [-r,T  + e].  An  application  of  Gronwall's 


■ 


inequality  immediately  gives 

Hxt  - x(e  J " c < L (6  ) h (c  ;T)etL(6 ) 
for  t€ [ 0 , T ] . This  shows 

lim  IIT^x  - Tt  x II  = 0 , x£  € , (3.11) 

e-»0  + 

uniformly  on  compact  intervals. 

Now,  choose  'xGX  and  a sequence  {x^l  in  If  such  that  x’^-»  “x  in 

X. 

Take  a sequence  e n -*  0. 

Assumption  (H5)  and  Proposition  1.2,  part  a),  together  with 
Remark  1.8  imply  that  for  any  u > 0 we  can  choose  a number  kQ 
such  that 

|x<e")(t)  - x(en)(t;xk)|  < y, 

|x(t;xk)  - x ( t ) | < y 

for  all  n and  all  t£[0,T]  provided  that  k > kQ.  Therefore  we 
obtain 

tx(en)(t)  - x ( t ) I < 2 y +lx(en)(t;xk)  - x(t;xR)l  for 
all  n, te[0,T]  and  k > k . 

By  (3.11)  for  each  k we  can  choose  nk  such  that 
1 x nk\t;xk)  - x(t;xk)  I 5 -jL 

for  telO.T] . Therefore 

, ( enu ) ^ 

limx  (t;x)=x(t) 

k-»o 

uniformly  on  compact  intervals.  Since  en  -»  0 was  arbitrary  this 
proves  the  desired  result. 
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Corol 1 ary . Under  the  conditionsof  Proposition  3.3  Ttx,x£X, 
is  the  unique  solution  of 
t 

Ttx  = Stx  + JSt_s  Aj  T$xds  , t > 0.  (3.12) 

Proof . In  view  of  formula  (3.8)  and  the  uniform  convergence 
of  T (€  ^ x on  compact  intervals  we  only  have  to  prove 
t , v t 

lim  JS,  A,  T^eJx  ds  = / S.  CA,T  x ds.  (3.13) 

e->0  + o 1 ,e  s o t"s 

First  we  consider  x£  C . Then 

IIAl,c  TiC)>r  - A1TS7"  ' lf£<Ts°*>  ■ f(Tsx)l 

■ lf<s,(s)Ti°*>  - 

< L(«)IISt(s)T(')x  - Tsxlle 

s L(6>{»s,(s)Tse)y  - s,(s)Ts^"e  + lls,(s)Ts?  - Ts*«e 

where  t(s)£[0,e]  and  6 = sup  ( NT  ^ e ^x  IU,  IIT  x IL)  . 

[0,T]  s esc 

This  estimate  shows  that 

lim  ||A,  T<e)x  - A,T  x||  = 0 
e-*0  1»e  s 1 s 

uniformly  on  [0,t).  This  proves  (3.13)  in  case  x£  if 

Now  take  x£X  and  choose  a sequence  ( x k ) in  if  such  that  xk  -»  x. 

For  each  k we  have 

t 

Tt*k  = stV  /st-sAiTs*kds- 

Recalling  the  definitions  of  and  we  see  that  for  y£X 

st-sAiTsy  = (f(Tsy>’  *s(y)) 

where  ips(y)  denotes  the  function  in  Lp ( - r , 0 ;R  ) defined  by 


*s(y)(0) 


f f(Tsx) 
( ° 


for  O > s-t , 
for  0 > s-t . 


Therefore 

/<St-sAlTs*  - St-sAl<s*lc>ds 

t t 

= (/(f(Tsx)  - f(Tsxk))ds,  J(*S(X)  - ^s(xk))ds). 

For  the  first  component  we  get 

t t 

l/tf(T  x)  - f(Tsx,)]ds|  < YB(t)(/lx(s;x)  - x(s;x.)lPds) 
o p - r 

where  Y^(t)  is  appropriately  chosen  according  to  (H2). 

For  the  second  component  we  have 
t 

ll/Us(x)  - <PS  (xk)  Ids  llp 

< /(  /(*s(x)(0)  - *s(xk)(0)  I pd0  ) 1/p  ds 

o -r  s s K 

.<  /(  /Mr(Tsx)  - f(Tsxk)lpd0)Vp  ds 
o -r  s s K 

< r^p  Yp(t)(  Jlx(s;x)  - x ( s ;xk ) I pds ) ^p 

The  estimates  (3.14)  and  (3.15)  together  with  Proposition 
1.2  prove  that 

t t 

Jim  /VsAlTs*kds  = / St-sAlTs"ds- 
k-*»  o o 


This  finishes  the  proof  of  the  corollary. 


3o  - 


Ik 


■ 


i 


& 5.  Representationof  the  solution  semigroup,  local  case 

In  this  section  we  remove  the  assumption  that  (H5)  is  a global 
condition  and  show  that  the  results  of  Section  3 remain  valid 
with  obvious  modifications. 

The  operators  A , A,  and  the  map  f are  defined  as  in  Section  3. 

€ 1 9 c £ 

We  need 

Lemma  4,1 . For  each  3 > 0 and  p > 0 there  exists  a constant 
*0  = *0(3>p)  such  that  for  all  A€(0,Ao)  and  all  x'eX  with 
II x||  < 3 the  equation 

(I  - AAJy  = x (4.1) 

has  a unique  solution  y€<$  such  that  llyll  < 3+p . 

Proof . The  equation  is  equivalent  to  the  fixed  point  equation 

y = ( I-XAo)_1x  + X(I-XA0)-1A1>ey. 

Then  the  result  follows  by  standard  arguments  using  the  dissi- 
pativeness of  Aq  and  the  Lipschitz  continuity  of  A^  £. 

If  for  x'eX  there  exists  a ball  with  center  in  the  origin  such 
that  equation  (4.1)  has  a unique  solution  y in  this  ball  then 
we  put  y = (I  - xA  )-1x.  The  lemma  shows  that  (I  - xA  )"*x 
exists  for  x sufficiently  small.  In  general  equation  (4.1)  does 
not  have  a unique  solution  in  all  of  3. 

Proposition  4.1.  Suppose  that  (HI)  - (H3)  are  true  and  define 
T | e ^ x”  = ( x ^ e ^ ( t -,x) , x|e  ^ ( x ) ) where  x^^tjx")  is  the  solution 
of  (3.6)  with  initial  data  x”.  Then 

t[e)x  = lim  (I  - £Aefnx  (4.2) 

n-*oo 

for  all  x'eX  and  te[0,t~)  the  convergence  being  uniform  on 

Jm 
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n is  globally  Lipschitzean  with  Lipschitz  constant  2 (cf.  [ 1 0 3 ) - 
Therefore  is  globally  Lipschitzean  with  constant  ( e ) . 

If  we  define  ft,  x = (f  (jT),0)  and  = A + ft,  , then 
ft  - al  is  dissipative  with  a = - + By  Remark  3.4 

C P o p 

we  see  that  Proposition  3.2  applies  and  we  get 

t[°x  = lim  (I  - )"nx  (4.4) 

n-*°o 

uniformly  on  [ 0 , T ] . Note,  that  the  solution  x^^  (t’,x)  of 
x(t)  = f ( x ( t ) , x t ) coincides  with  x ^ e ^ ( t ; x ) on  [ 0 , T ] . 

In  order  to  finish  the  proof  we  have  to  show  that  in  (4.4)  we 
can  replace  (I  - ^e)~n  by  (I  " ^Ae)~n  provided  that  n is 
sufficiently  large.  To  this  end  we  only  need  to  prove  that 
11(1  - ^e)"^x||  < 3 , j=l,...,n,  for  all  te[0,T]  and  n > nQ, 
nQ  some  positive  integer  larger  than  the  constant  Xq(3,1)  of 
Lemma  4.1. 
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We  first  suppose  that  xe^)  and  write 

0 - - T<£)  x - 


M - Jift  \ ~ J +K  ( j ( e ) x'  - n - l.ft  \ T ( e ) 
nMc  ' UKt/n  x nV  1 (k  + 


i (i  - j - (I 

K =0 

■ V(i  - |A,)‘J+,ctjn(A  T<‘> 


(k  +1 ) t/n 


ET(Kil)VnX  - AeTKt/n  +sx>ds 


This  and  the  dissipativeness  of  - al  give  the  estimate 

"<<  - S '“‘j’  tf“*eT(‘Il)Vn*-*tTiwn*ll<ls 

From  Proposition  2.3,  part  b)  there  exists  for  each  positive 
number  p a number  nQ  such  that 

lift  J ( E ) % _ ft  T ( E ) x’ll  < ii 

e 1 (K  + l )t/n  X e Kt/n+s  x 11  < u 
for  n > nQ  , se[0,t/n]  and  k =0 , . . . , n- 1 . Thus  we  arrive  at 

||(I  - ^e)“Jx  - T j ^nx  II  < u eat  < pteat  < 1 

for  n > n , j = l,...,n  and  t£[0,T]  provided  that  p is  chosen 

- 0 t 

such  that  pT  e°  < 1.  By  choice  of  3 the  last  estimate  proves 
||(I  - dx||  < 3 on  [0,T]  for  n > nQ  and  j = l,...,n. 

If  x is  not  in  the  result  follows  from  the  densitiy  of  3 in  X 
and  the  Lipschitz  continuity  of  (I  - 

Proposition  4.2.  Suppose  that  (HI)  - (H5)  are  true.  Then  for  all 


1 im  T ^ e ^ x*  = T . x 
e-»0+  z z 


uniform  on  compact  subintervals  of  [0,t~).  Moreover,  T^x  is 

the  unique  solution  of 

t . 

Ttx  = Stx  + Jst_sAiTs*  ds  , t+[0,t~). 
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Proof . We  firs1-  prove  that  for  T£(0,t~)  there  exists  a number 
cQ  > 0,  such  that  the  solution  x^c^(t;x)  of  equation  (3.6)  with 
initial  data  7 exists  on  [0,T]  for  0 < e < e . Moreover,  the 

( p \ “0 

functionsx'  '(t;x)  are  uniformly  bounded  on  [0,T]  for 
0 < e < eQ.  Indeed,  choose  6 such  that  IIT^.x'll  <6  - 3/2  on 
[ 0 , T ] . Assumption  (H5)  and  Proposition  1.1  imply  that  there 
exists  a constant  a > 0 such  that  x^e^(t;x)  exists  on  [0,a] 


exists  a constant  a > 0 such  that  x^e^ 


and  ||T  £ 'x^l  <6-1/2  on  [0,a]  for  all  e > 0.  Note,  that 
f (0,0)  = f(0,0).  Then  quite  analogous  to  the  proof  of 
Proposition  3.3  we  see  that  -*  T^x  uniformly  on  [0,a]. 

Therefore  there  exists  a number  > 0 such  that  IIt|e^x||  <6-1 
on  [0,a]  for  0 < e < e^.  Again  by  (H5)  and  Proposition  1.1  we 
obtain  that  x ^ e ^ ( t ;t!.c  ^ x ) = x^e^(t+a;x)  exists  on  [0,al  for 


0 < e < £.,  and  ||T v 'x 
- 1 a + t 


x)  = y(t+a;x)  exists  on  [0,a]  for 
<6-1/2  on  [ 0 , a ] . Then  we  can  chose 


a number  62  < 


such  that  ||T. 


<6  - 1 on  [0,2a]  for 


0 < c < c 2>  In  a finite  number  of  steps  we  arrive  at  a number 
> 0 such  that  T^x  exists  on  [0,T]  with  ||t{c^>T||  <6-  1 on 


[0,T]  for  all  t € ( 0 , e ] . The  rest  of  the  proof  is  quite  ana- 
logous to  that  of  Proposition  3.3. 

Remark  4.1.  If  [0,t~  _)  denotes  the  maximal  interval  of 
( e 1 ~e 

existence  of  xv  '(tjx')  then  the  proof  of  Proposition  4.2  shows 
that 

1 imi nf  t£  > t~. 
o _ rw  x,e  - x 
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& 6.  Numerical  approximation  by  averaging  projections 

The  objective  of  this  section  is  to  study  the  approximation  of 
solutions  of  the  nonlinear  f u nc t i ona 1 - d i f f eren t i a 1 equation 
(2.1)  by  solutions  of  ordinary  differential  equations  acting 
on  a sequence  of  finite  dimensional  subspaces  of  X.  We  take 
the  averaging  projections  of  Banks  and  Burns  [2  ] and  deal 
with  the  case  where  (2.1)  is  the  difference-differential 
equation 

x(t)  = h ( x ( t ) ,x(t-k1) , . . . »x(t-km) ) , t > 0,  (5.1) 

0 = ko  < kl  <•••  < kpi'  assumptions  on  h are  the  same 

as  stated  in  Remark  1.2.  Equation  (5.1)  is  the  more  difficult 
case  compared  with  an  equation  where  a condition  of  type 
(1.5)  holds.  Some  technical  difficulties  which  we  couldn't 
overcome  are  the  reason  that  we  finally  can  get  the  desired 
approximation  results  only  for  initial  data  x = (cp(G),(p), 
where  ip  is  Lipschitz  on  [-r,0],  i.e.  <p€W  * ’“(  - r , 0 ; IRn  ) . 

We  denote  the  subspace  of  these  initial  data  in  X with 
‘ftf’  *°°.  For  the  rest  of  this  section  we  fix  a solution 
x(t;x),  xe?<yl,0°’  of  (5.1)  which  exists  on  [-r,T],  T > 0. 

There  exists  a constant  K such  that  I x ( t ; x")  I < K on[-r,T]. 

We  define  the  map  h(yo>...,ym)  = h (yQ  , . . . ,ym  ) where 

Yj  = yj  for  I y j I < K+p  and  = yy^y  y.  for  I y j I > K+p . 

J 

The  positive  number  p will  be  chosen  appropriately  in  the  course 
of  the  proof.  It  is  easily  seen  that  ft  is  globally  Lipschitz 
with  Lipschitz  constant  ft  = L(K+p).  According  to  Remark  1.2 
we  can  take  X = (RnxL^(-r,0;Rn)  as  state  space  for  the 
equation 


x(t)  = ft(x(t)  ,x(t-kj ) , . . . ,x(t-km) ) 

By  definition  of  ft  it  is  clear  that  x(t;x)  is  also  the 


(5.2) 
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solution  of  (5.2)  with  initial  data  x on  [-r,T]  . 


N 

The  sequence  (X  ) of  subspaces  of  X is  defined  as  in  [2  ], 

N N n N 

i.e.  a typical  element  of  X is  ( n , I v • x^)>  where 

J J 


j = l 


N N n J N 

n»Vp...,V|y  are  elements  in  R and  Xj  is  the  characteristic 

function  of  the  interval  [t^,t^_,),  j=l,...,N  , t^  = -j£. 

J J J-  J " » 


N 


,N 


N, 


N 


N N. 


The  projections  n : X -»  X are  defined  by  it  (n,<p)  = (n  » I ip^x,) 

j = l J J 

where 


tN 

= J /3“1*p( s)ds  , j = 1 N . 

J tj 

N , 


(5.3) 


It  is  easy  to  see  that  tt  is  a contraction  on  X for  N = l,2 

Moreover , 


limn  x = x , x£  X 
N-*» 


(5.4) 


If  we  define  ll(n»tp)ll0O  = max  ( I n I , ll<P  11^)  for  (n,<p)€R  *L  (-r,0;R  ) 
then  an  easy  calculation  shows 


llnN(n,<P)  11^  5 H(n  »ip)  11^ 

for  (n,<p)€lRnxLoo(-r,0;Rn)  . Note,  that  IIttN  ( n ,cp)  II  = 

N N N °° 

max  ( | n | , I (pj  I , . . . , I I ) , where  the  cp  j are  defined  in  (5.3). 


Following[2  ] we  define  the  operators  : X 


N . N 


N 


XN  by 


■yN , ^ N N . , n N ^ , N , 

= ° = l VJ  “ 1 ‘ 


where  v^  = n and  put 


A^x  = S’q  tt^x  , x£X 


(5.5) 


(5.6) 


(5.7) 


N < N 

Obviously,  Aq  is  a bounded  linear  operator  with<</(A0)  = x 

and  therefore  generates  the  linear  CQ-semigroup 
sj  = eAot  , t > 0. 
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N N 

We  need  some  results  on  Aq  and  S^.  with  respect  to  the  space  X 


m i 

Lemma  5.1.  a)  The  operator  Aq— ^1  is  dissipative  for  all  N and 


lisjxll  < e1/2  t||x||  , x6X  , t > 0. 


b)  For  all  x6X 


1 n IN  ~ r'  ~ 

1 im  Sf x = Sf x 
N-« 


uniformly  on  compact  intervals 


Proof . In  order  to  prove  part  a)  we  f i r-.t  take 

x = (n,  Z From  (3.3)  we  get  (note,  that  the  duality 

j = l J J 

map  is  the  identity  map  in  case  p=2) 

~v  _ " , N . N . N 


✓ . IN~  ~V  /IN  IN  . IN 

<A0x,x>  = I (Vj.j  - »j)Vj  = 

J A 

1 N N N 2 1 2 

= ' Vj'  + 7lnl 

1 2 1 ~ 2 
< £lnl  < £llxir  • 


1 . N . 2 
? VN 


For  general  x”  = (ri,cp)eX  we  write  x = n^x'+y^  where  n^y  = 0 and 
therefore  A^y^  = 0 according  to  (5.7).  This  gives 


<^A^x,x)>=  ( A^ttNx  , ttNx  ) < I n I 2 f-^llxll' 


N N 1 

The  estimate  on  St  is  clear  by  the  dissipativeness  of  A --^I. 

Part  b)  is  clear  by  standard  results  of  linear  semigroup  theory 
(cf.  for  instance  [15]).  We  refer  to  [2  ] for  a proof  of 
A^x  -*  A„x”  for  xe<0  . 


Remark  5.1.  If  we  choose  the  columns  of  the  matrix 

( I , I ( s ) , . . . , I x!!  ( s ) ) as  a base  of  XN  then  has  the  following 
I in  n o 

matrix  representation  [5^]: 


. N r-vN  , N 

w = [A  ]w 

with  w^(0)  = (n  » v j , . . . , vjj ) . For  x = nNx+yN  we  get  S^x  = 

S^n^x’  + S^y^.  where  S^y^  = y ^ . By  (5.4)  and  part  b)  of  Lemma  5.1 
we  get  yN  -»  0 and 

lim  S^tt^x’  = , x£X , 

t-*oo 

uniformly  on  compact  intervals. 


Remark  5.2 . Since  [ ^ ] has  a special  structure  it  is  not 
difficult  to  obtain  an  explicit  representation  of 
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We  need  the  explicit  representation  of  [e  ] in  order  to  prove 
N 

Lemma  5.2.  S^  , t > 0 , is  a uniformly  equ  i co n t i nuou s semigroup 
on  RnxL°°(-r,0;IRn),i.e. 

Iisjjxll^  > 311x11^  , x€lRnxLoo(-r,0;Illn)  , 

for  t > 0,  N = 1,2 

Proof.  Since  S^x”  = S^n^x  + y^;  by  Remark  5.1,  and  therefore 

IlS^xll  < IIS  ^ tt  N x II  + lly N II  < IIS  ^ n N x II  + llx  II  , we  only  need  to 

N N 

prove  that  St  is  uniformly  equ  i conti nuous  on  X with  respect 

L ~ N N N 

to  the  sup-norm.  For  x = (n».l  v^Xj)  we  get 
[a^o t]  col  (n.vj v{[)  = 

7TN  t rNt  N N 

[eMotj  col (n.0, . . . ,0)  + [eA°r ] col (0, v J, . . . , vj) . 

From  the  representation  of  exp[AQt]  given  in  Remark  5.2  we  get 

a = [eAot]col(n,0,...,0) 

Nt  N+  Ni  N-l  , N.  . 

= col  (n  » ( l-e~r  ( 1-*—) ) n »•••.(  1-e" r I |r(^)J)n) 

r j=0  r 


b = [e^t]Col(0,v5l,...,v[j) 

. . N_t  Nt  N-l  . ..  . . . . 

= col  (0,  v^e”r  ,...,e-^  Z^  yr^.  j (^) 3 ) • 

From  these  explicit  representations  we  immediately  get 


lla  ||  < I n I . lib  II  < max  | v • 

CO  - . J 


which  proves 


||[e^ot]col  (n.vj,  • • . , v JJ ) 11^  < 211x11^. 


For  x 


(n»<p)€X  we  define 


?N(x)  = h(n,  (PN(-k1) , . . . ,ipN(-km) ) 

a nd 

fN(x)  = h(n,4>N(-k1),...,cpN(-km)) 

N N ~ 

where  (n,<P  ) = n (n,4>).  Note,  that  and  fN  are  well- 
defined  on  X since  is  right-hand  continuous  on  [-r,0). 

The  corresponding  definition  in  the  general  case  of  equation 
(2.1)  is  f^  = fonN  which  makes  sense  if  XNc^(f)  for  all  N. 

As  an  immediatly  consequence  of  (5.5)  and  the  Lipschitz- 
continuity  of  h and  h we  get 

1 f M (y ) I 5 L(M)  Nx-yll^ 

for  “x  ,y£  Rnx  L°°(  - r , 0 ; Rn  ) with  llxll^  < M,  llyll^  < M and 

l?N(x)  “ V?)1  - ^llx-yl^  (5.9) 

for  x,yeiknxLoo(-r,0;IRn). 


Since  ||nNx  - nNy|loo  < (£)V2  llnNx-nNy||  < (^)V2  ||x-y||,  we 
also  get 

irN(x)-rN(y)l  < (^)VZ  Cllx-yll 


(5.10) 


for  x,yeX 


In  order  to  make  the  representation  not  to  complicated  we 
only  deal  with  in  detail  and  give  the  modifications  for  f^ 
in  a remark.  Corresponding  to  we  define  the  nonlinear  operators 


4o 


, CN/ 


A”x  = (fN(x),0)  , xex 

and  AN  = A^J+aIJ* . 

o 1 

N 1 N 1/  2 ** 

Lemma  5.3.  a)  A - .ol  is  dissipative  with  o>  = -£+(—)  *-  • 

N 

b)  A is  the  infinitesimal  generator  of  a strongly  continuous 
semigroup  , t > 0,  such  that 

lujx  - Tjyll  < eu)tllx-yll  , x,y€X  , t > 0. 

N~ 

Moreover,  T.x  is  the  unique  solution  of  the  integral  equation 


tN~  cN~  }cH  .NtN~.  . n 
T^x  = S^x  + JSt_sA^T^xds  , t > 0. 

o 


(5.11) 


Proof . Part  a)  follows  easily  from  Lemma  5.1,  a)  and  (5.10). 

Part  b)  is  proved  by  an  application  of  the  Crandal 1 -Liggett 
theorem  [9  ] or  follows  from  [18]  where  also  the  result  re- 
garding the  integral  equation  (5.11)  is  proved.  If  one  applies 
the  Crandal 1 -Liggett  theorem  one  has  to  prove 

Ti(  I - x A N ) = X for  X£(0,i).  But  this  is  easily  done  by  a discussion 
of  the  fixed  point  equation 


y - (i-xajfhr  * xd-xA^-^Jy  , x ,y€X 


N 

Remark  5.3.  We  choose  the  same  base  of  X as  in  Remark  5.1  and 
define 

FN(n,Vj, . . . ,vj)  = col(?N(n,  I vjXj)>  0,...,0). 

N J = * 

~ '’WWW  M 

For  x = (n>  I v-x-)€X  let  w (t)  denote  the  coordinate  vector 
j = l J 3 

of  T^x*.  Then  wM(t;x)  is  the  solution  of 


wN  • tsr"  iw"+  rN(w")  . t » 0 , 


rN,..N 


N. 


(5.12) 


N N N 

w (0)  = col(n,v  ,...,vw).  This  follows  from  the  integral 
equation  (5.11)  and  the  fact  that  X is  invariant  with 
respect  to  AM.  For  general  xgX  we  write  >T  = n^x'+y^  and  get 
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N~  N N~  ut  ~ N~ 
lITjx  " Ttn  xN  - e NX-"  xll  , t > 0, 

i.e.  by  (5.4)  we  have  lim  IIT^x  - T^ttNx||  = 0 
uniformly  on  compact  intervals. 

Corresponding  to  equation  (5.2)  for  e > 0 we  define  f as 

£ 

in  Section  3 and  get  from  Proposition  3.3 

lim  T^e^x'  = T.x  , x€X  (5.1 

e-»0+  1 1 

/ £ \ 

uniformly  on  compact  intervals,  where  T£  ' , t > 0 , is  the 
solution  semigroup  of  equation  (3.6)  with  infinitesimal  gene- 
rator A = A + A,  , A,  x = (f  (x)  ,0)  . The  calculations 

given  in  Remark  1.2  show  that  the  function  y (t)  in  (H2)  can 

~ 1/2 

be  chosen  independently  of  3 , y ( t ) = (m+l)L  t . This 
function  can  also  be  used  for  all  equation  (3.6)  wi  th  e£(0 , ki  ] . 
Therefore  Af  - ol  is  dissipative  with  (->=2+jy(e)  = -^+ (m  + 1 )Te" 
by  Proposition  3.2,  e < min(kj,l). 

In  the  next  step  we  define  for  e > 0 


fN,c^X^  " e /VSxX>dT  ’ x6X> 

and  the  operator  A^  = A^  + A^  ^ where  A^  £x*  = (f^  e(x'),0)  for 
x e X . 

Lemma  5.4.  a)  A^  - ul  is  dissipative  with  o = + (m+1) 

for  all  N > N and  ANx*  -»  A x'  for  x€ 
o z z 

N 

b)  Ae  is  the  infinitesimal  generator  of  a nonlinear  semigroup 


(e,N) 


, t > 0,  such  that 


(e,N)~  x(e,N)~..  ut 

t x “ T t yl1  - 6 »x-yii 


for  t > 0 x ,y£X  and  N _ N . 


c)  T £ e ’ ^ x is  the  unique  solution  of  the  integral  equation 


tU,N) 
1 1 


;N  ,N  t(c,N) 


* ■ * / s,;-s  svr 


xds  , t > 0 


(5.14) 
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d)  For  all  e > 0 we  have 

lim  t{G,N)x  = tJe)x  , x£X, 
uniformly  on  compact  intervals. 

Proof.  We  first  prove  that  fw  is  globally  Lipschitz  on  X with 

E I/O  fci 

Lipschitz  constant  (m+l)te  , which  proves  that  A - wl 


Lipschitz  constant 
is  dissipative. 


For  x.j  = (o^,(P^)GX  , i = 1,2,  define  the  function 


(s) 


i ( s ) on  [ - r , 0 ) , 
ni  on  [0,”)  , 


, i = 1,2  . 


Then 


1 ^N,e  (X1  )"fN;e(x2  ) 1 - (n  1 ’ 1 )"^N  ^ 2 ’ ^ 2 ^x  ^ dT 

tN  _ ° tN 

= 1 <P!  ( s+t  )ds, . . . ,£j  m tp1(s+x)ds)  - 


- h(n2 - 


J 1 


(P2(s  + t )ds)|  dx 


where  we  choose  N such  that  J < min(k.  k.)  and  j is  such  that 

NN  ' j J + J ^ 

k^Gttj  ,tj^_j).  Since  TT  is  globally  Lipschitz  we  get  (for  e<  1) 

~ r m 


|fN,e(xl)  ' f| 


r f f N m P"  \ 

I , e ( X2  > 1 i "e  J\lnl“n2l  +rI1Jw  I ( ^ + s ) -<p  ( T + s ) I ds). 

0 »M  tH 


tN 

l.J  -1  e 


- ?(?  J I ^ i “ ^ 2 I c*sdl  + 7 1 f J Ii1'1(t  + s)-^„(t  + s)  Idsdi) 

\ tN  o u=i  *;n  o £ / 


I t 
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<f^eV2(/  ( }l^1(s)-^2(s)l  2ds ) 1/2  di  + I ( /l*1(s)-*2(s)  I 2ds)V2ch; 

e r tj  -r  1 £ y = 1 tj  -r  1 * 


m r y 


t" 

J -i  0 


< te"V2  (m+l)||xrx2H. 


In  order  to  prove  the  second  claim  of  part  a)  for  7=  (<p(  0 ) ,ip)€cv  we 
consider  the  estimate 

e 


|fN>c(x)-fe(x) ! < ~ /l^N(STx)-f (Stx)  Idx 


t? 


N Jj  1 


-1 


= ~ / I n (tp(  0)  ,-p  j N cp(t+s)ds,  . . . )-h(ip(0),(p(x-k1),.  . . )]  di 


t'; 


tN 

t r N.  - A’1 


1 


JN  l«p(x+s)-«p(T-k  )|ds)d 
o y = 1 t, 

Jy 


where  we  again  choose  N such  that 


H < ">1n(kj  + 1-kj) 

J 


Since  <P  is  uniformly  continuous  on  [ - r , 0 ] for  arbitrary  y > 0 


we  can  define  N such  that  for  N > Nn  we  have 
o - o 


N J 


I <p(  t + s ) - tp(  t-  k )l  < y for  sG[ti,t_1  ]. 

P J J i 


From  the  estimate  given  above  we  therefore  obtain  (note,  that 

k eltj  , t ^ . ) ) for  N > N 

y J J - 1 ' ' 0 

y y 

1 - 7 7"R  p£  = 


which  proves  A,  x -♦  A,  xfor  xGc3  as  N 

lit  I , e 


The  proof  of  parts  b)  and  c)  is  quite  analogous  to  the  proof 
given  for  Lemma  5.3. 

Part  d)  follows  from  Corollary  4.2  of  [8  ] using  part  a)  of  the  lemma. 
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Lemma  5.5.  For  each  x€u/*  ’°°  we  have 


lim  ||T 
e->0  + 


<£'N>x  - All  - 0 

L L °° 


uniformly  with  respect  to  N and  uniformly  on  compact  intervals. 


Proof . From  (5.11)  and  (5.14)  we  get 

IIT<e’N>x  - TNtxlL  < /IIS?.s(fN,£(T<E,N>x)-r  n<T«x)  ,0)  ll^ds 

A * • 


5 3/'fN,e(Tse,N)^)  ' Tn(T^x) ! ds 
o 


f (TN; 
N S 


5 3/lfN>c(T<£'N)x)  - fN>£(T|!x)lds 
+ 3^|fN,e(Ts^)  “ ^(T^)lds 


N v s' 


< 3r;iu(e’N)x  - 1H  Xll^ds  + 3/ 1 f N (tJx)  - fN(Tjx)|ds.  (5.15) 
0 0 ’ 


Here  we  have  used  the  fact  that  from  (5.9)  we  immediately  get  the 
same  condition  for  all  f^  . For  the  second  integral  on  the  right- 
hand  side  of  the  estimate  we  get 


/'Ve(Ts*>  ' rN<TsX>'ds 

< 7 J Jll»NStT^x  - tJx|L  drds. 


(5.16) 


N N~  ^ N N 

Using  for  fixed  t the  representations  n S T x = (n(s),  I ij> . ( s ' x •; ) 

N~  N m n T s i=l  J j 

and  T%  = (n(s),  X «pj(s)x?)  we  get  J 

j = l J J 

||nN  St  T^x  - TgXll^  = maxltPj(s)  - j ( s ) I . 

J 

Recalling  the  definition  of  we  get 

, N N N,  . N , . N . w N N , 

♦j  - 

N / . N j . N \ / N N, 

* ?(tj-l+'  - tj-x>ld>j.1<-  Vj) 


(5.17) 
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N N N 

where  < is  defined  by  t-+T£[t.  , , t - .)  and  we  have  dropped 

J J-K  + l J " 

the  argument  s. 

For  any  natural  number  we  have 

tN 

liPj_p-  Vjl  < 7 /Jj’1l‘P(s  + pjf)-<p(s)  Ids 
tN  s+pr 

< 7 /^_1  J " fT I 4> ( o ) I d a ds  < p5,*,L 

V 

Therefore  we  get  from  (5.17) 


5 [(tj-K_(tj+T))(K"1)  + (tj-l+T_tj-K)K]l|(i)|1co  = 


This  estimate  together  with  (5.16)  shows 

£ lfN>C(Ts*)  ‘ fN(Ts^|ds  ^ ZztWvll 
Using  this  in  (5.15)  an  application  fo  Gronwal's  inequality  gives 
l)T(e,N)~  _ tN~|^  < 3^e  t ||(p||ooexp(3rt) 
which  proves  the  Lemma. 


We  now  are  in  the  position  to  prove  the  main  result  of  this  section 
Define  corresponding  to  f^  as  F^  was  defined  in  Remark  5.3 
correspond i ng  to  TN . 

Proposition  5.1.  For  ’°°  let  u^(t)  be  the  solution  of 


Then 

lira  xN(t)  = Tt x"  . 

N-k» 

uniformly  on  compact  intervals. 

Proof . We  first  deal  with  the  global  case,  i.e.  with  equations 
(5.2)  and  (5.12). 

For  y > 0 according  to  Proposition  3.3  and  Lemma  5.5  there 
exists  an  eQ  > 0 such  that 

HTtx  - t[£o>xII  < y , 

HI  ( eo  » N )~  _ TN~n  < w 

for  t€  [ 0 , T ] and  all  N.  Then  we  can  choose  Nq  such  that  for  N > Nq 
and  t G [ 0 , T ] (cf.  Lemma  5.4,  d)) 

HTjox  - t[co»n)x||  < y 

The  estimates  together  give 
l!Ttx  - Tjxll  < 3 y 

* for  all  N > N and  t£[0,T).  This  and  with  Remark  5.3  proves  the 

result  for  the  global  case. 

If  we  consider  the  local  case  we  choose  the  constant  P which 
appears  in  the  definition  of  ft  equal  to  1 and  get 

N 

Ttx  = li m (wj[(t),  Z ( t ) x?)  (5.19) 

1 N-°o  0 j = l J J 

M 

uniformly  on  compact  intervals,  where  w (t)  is  the  solution 
of  equation  (5.12)  with 

N 

N w N N N N~ 

wN(0)  = col  (g)(0)  ,<Pj  ....  ,ipN)  , (<P(0),  I <PjXj)  = n x- 

J — * 

we  only  need  to  show  that  for  N sufficiently  large  we  have 

j 
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HwN(t)lloo  < K+l  on  [0,T].  Then  w^(t)  = u^(t)  for  N sufficiently 
large,  because  ?^(u)  = F^(u)  for  llull^  < K+l. 

Using  Lemma  5.2  and  (5.9)  we  get  from 

wN(t)  = [e^ot]wN(0)  + / [e^o(t_s)]  T (wN(s))ds  (5.20) 

o 

by  an  application  of  Gronwal's  inequality  that  there 
exists  a constant  P 

llwN(t)lloo  < P 
on  [0,T]  for  all  N. 


Differentiating  (5.20)  we  get 

.*Xt 


wN(t)  = [eAot][^]  wN(0)  + /[eAo(t_s)](^]FN(wN(s))ds 


+ FN(wN(t)) 


N 


(5.21) 


An  easy  calculation  shows  that 

ll[^]wN(0)li  = IIaJxII  < P,  for  all  N 
where  P^  = llipll^  (note,  x = (ip(0),tp),  ipe/tT  ’°°)  . 
Therefore  according  to  Lemma  5.2 
Ille^nffJjj  wN(0)lloo  < 3Pj 


If  we  write  the  second  term  on  the  right-hand  side  of  (5.21) 
explicitly  we  get 


7 /col(0,w-"(t-s)r„(wN(s)), 


1 


• ’(Y-Typ  r 


(Xt^sJN)N-le-r(t-s)r  (wN(s)))ds  > ,(5<22) 


Using  the  uniform  Lipschitz  condition  on  and  the  boundednes 

N M 

of  llw  (t)ll  we  after  some  elementary  calculation  see  that 

there  exists  a constant  P,.  such  that  the  sup-norms  of  (5.22)  is 

uniformly  bounded  by  P^.  It  is  clear  that  the  third  term  on  the 
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T 


right-hand  side  of  (5.21)  is  also  uniformly  bounded.  Therefore, 
there  exist  a constant  Q such  that 

llwN(t)llco  < Q 


on  [0,T]  for  all  N.  From  this  we  get 

1 wj_i ( t)  - wj(t),  = Jl Wj (t) I < JQ 

and 

l"j+i<t)  - 1 = ^ 

for  all  N and  t€ [ 0 , T] . 

For  any  natural  number  p we  therefore  have 

|Wj-P(t)_WJN(t)l  - plJQ’ 

|WjN+p(t)  - WjN(t)|<  p-^Q 


(5.23) 


for  all  N and  t£ [ 0 ,T]  . 

Suppose  that  for  t€[0,T]  there  exists  a subsequence  (N^)  which 
for  simplicity  we  again  denote  by  N such  that 

|w"  (t) I > K+l  (5.24) 

JN 


for  at  least  one  » 1 < < N.  By  (5.23)  and  (5.24) 

we  get 


|w"  (t)|  > K+* 


(5.25) 


N-p 


N 


for  all  k and  all  p < (to]  denotes  the  largest  integer 

<a).  If  Nk  denotes  the  smallest  integer  such  that  k = 
hen  < kN^  and  k ] > k.  Therefore  from  (5.25)  we 


then 
infer  that 
N 


ik  w^(t)x-k(s)|  > K+i  , k = 1,2, 
j = l J J 


(5.26) 


for  s in  an  interval!  of  length  greater  than 


V V 


r - 49  - ® 

t I 

I In  this  case  (5.26)  shows  that  we  cannot  have 

II  x,(t;x)ll2  - 0 

j=l  J J 1 

as  k -*  This  contradicts  (5.19)  and  therefore  llw^(t)lloo  < K+l 
for  all  N and  t€[0,T]  . 

Remark  5.4 . The  proof  given  for  Proposition  3.1  in  the  local 
case  can  easily  be  adapted  to  show  that  for  xek'  ,ot>  we  in  fact 
have 

limllx^(t)  - T.xll  =0 
, N-oo  °° 

uniformly  on  compact  intervals. 

Remark  5.5.  It  is  clear  that  the  results  of  this  paper  cover 
also  the  case  f^cr^f  )cMxL^( -r  ,0  ;M)  where  M is  an  open  subset  of 
Rn  and  fff/f  = { ( ip  ( 0 ) , cp ) | cp  6 C(-r,0;M)}. 

; Remark  5.6.  The  approximation  scheme  given  in  Proposition  5.1 

has  been  tested  for  various  nonlinear  examples.  In  each  case  the 
numerical  results  indicate  linear  convergence  uQ(t)  -*  x(t;x) 
as  N -*  oo. 
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Summary. 


In  this  paper  we  investigate  the  ontimal  control  problem 
for  Volterra  integrodif ferential  equations 


x(t)  = A (t)  x (t)  + |f  (t-s)  x (s)  ds  + B ( t)  u (t)  , 


where  the  target  sets  are  elements  of  some  function  space. 


: (t)  = i^(t)  for  0 < T-r  ^ t ^ T . 


The  approach  used  is  the  abstract  theory  of  Dubovitskii- 
Milyutin,  the  result  is  a necessary  condition  in  form  of 
a maximum  principle. 
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1.  Introduction. 


blern, 

*aJjc 


*In  this  paper  we  invcstfga-te^the  optimal  control 

gSU2_s  Xw/W  aaJjC- 

for  integroclif ferential  equations  of  Volterra  type^  where  the 
target  sets  are  elements  of  some  function  space.  The  approach 
used  is  the  abstract  theory  of  Dubovitskii  and  Milyutin,  the 
result  is  a necessary  condition  in  form  of  a maximum  principle. 

J V 

The  form  of  the  considered  equation  was  motivated  by  the 
applications.  A well  known  example  for  the  occurance  of  integro- 
differential  equations  of  Volterra  type  is  the  point  kinetics 
model  for  nuclear  reactors,  given  by  the  equations  (see  Akcasu 
et  al.  [1],  DiPasquantonio  and  Kappel  j.21]  ) : 

p(t)  =-4J(r)lP(t)  - VU1  " T2,$ 
r1  1 

(t)  = Aj[p(t)  - q^ft)]  , j = 1,2,...  ,N, 

where  p is  the  deviation  of  the  reactor  power  from  the  equi- 
librium level  n2  , q^  (t)  gives  the  j-th  group  delayed  neutron 
precursor  density,  1,  , Aj  , are  positive  constants  and  the 

reactivity  ^ = £(pfc)  is  feedback-functional  which  may  be 
taken  as 

i 

r 

£(pt)  = h*p (t)  + j f (t-s)p(s)ds  + u ( t ) , 

- eo 

h being  some  real  constant. 

Setting  x :=  col (p,q^ , . . . ,qN)  and  defining  A,  F,  B,  appropriate, 
this  system  can  be  written  as 

t 

r 

x(t)  = Ax(t)  + j F(t-s)x(s)ds  + Bu(t)  . 

-po 

An  important  problem  is  to  choose  u(t)  such  that  the  reactor 
is  shut  down  with  minimal  Xenon^S  build-up.  The  shut-down  state 
corresponds  to  x = col ( 1 , 1 , . . . , 1 ) . In  order  to  shut  down  the 
reactor  we  have  to  hold  x on  this  value  for  some  time  interval. 
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Another  problem  is  to  choose  u(t)  such  that  the  level  of  the 
power  equilibrium  state  is  changed  to  another  value  in  minimal 
time.  In  all  these  cases  a certain  terminal  function  has  to  be 
reached,  while  a certain  cost  functional  has  to  take  a minimal 
value. 

There  are  some  other  fields  where  Volterra-integrodif ferential 
equations  occur,  for  instance  viscoelasticity,  theory  of  material 
with  fading  memory,  population  dynamics,  spread  of  infectious 
deseases,  prey  predator  problems,  and  control  of  chemical  re- 
actions. 

\ 

| | 

We  shall  consider  the  equation 

t 

x ( t ) = A (t ) x (t ) + J F (t-  t)  x ( r)  di  + B ( t)  u (t)  (1) 

~OJ 

with  function  space  terminal  conditions.  This  equation  is  suf- 
ficient general  to  cover  the  above  mentioned  applications. 

In  chapter  2.  we  formulate  the  problem  and  the  main  result  of 
the  paper.  This  is  a maximum  principle  for  the  solution 
(xq,uo)  of  (1),  which  minimizes  a certain  cost  functional. 

In  chapter  3.  we  sketch  previous  approaches  to  problems  of 
such  kind,  which  gave  similar  results  by  using  other  methods 
of  proof.  Chapteis  4 to  1 1 present  auxiliary  material,  which  is 
needed  in  the  proof  of  our  theorem.  We  consider  the  "adjoint" 
equation  in  the  space  V of  functions  of  bounded  variation  and 
prove  the  existence  of  an  solution  to  this  equation.  In  chapter 
7.  the  theorem  of  Dubovitskii  and  Milyutin  is  cited  in  the  used 
form.  In  chapter  8.  we  discuss  the  spaces  that  may  be  used  and 
motivate  our  choice  of  the  space  C*5  (-<*>, ({]  of  continuous  bounded 
functions.  In  chapter  9.  we  give  an  existence  proof  for  solutions 
to  (1)  in  the  space  Cb.  Now  we  proceed  to  the  application  of 
the  theorem  of  Dubovitskii  and  Milyutin.  In  chapter  1o.  the  cone 
of  tangent  directions  to  the  solution  set  of  (1)  is  given.  Then 
in  chapter  11.  we  determine  the  dual  cone  K to  this  tangent 

cone.  Now  the  Euler  equation  of  the  problem  can  be  established 
and  is  used  in  chapter  12.  to  derive  the  desired  maximum  condition. 
We  conclude  with  some  remarks  concerning  further  necessary  work. 


L . „ - _ A 


2.  The  main  result 


Our  aim  is  to  prove  the  following  necessary  condition  for  the 
solution  of  the  below  defined  control  problem  (CP) : 


Theorem.  Let  xQ(t),  uQ(t)  be  a solution  of  the  control 
problem 


(CP)  : 


t 

x(t)  = A ( t ) x ( t ) + jF  (t- 1)  x (T ) dt  + B ( t ) u ( t ) , 

— CO 

x(t)  - ^(t)  for  -co  < t £ 0 and  T-r  & t 4 T , 
where  T-r  > 0 , 


l 

| <^(x(t)  » u ( t ) , t ) dt 


min  , 


where 

1 . A,  B are  nxn  resp.  nxr  matrices,  continuous  and 
of  bounded  variation, 


2.  F (s ) is  of  bounded  variation, 
o 

J |F  (s ) | ds  < oo  , 

4-  «> 


3.  is  continuous  and  bounded  on  the  intervals  -oo  < t ^ 0 , 
T-r  £ t & T. 

Now  let 

4.  '^(x,u,t)  be  continuous  in  x,u,  measurable  in  t, 

continuously  differentiable  with  respect  to  x,u,  and 

$ (x,u,t),  ^ (x,u,t)  be  bounded  for  bounded  x,u, 

X u 

r 0 

5.  M be  a closed  convex  set  in  R , with  M £ 0 , u(t)  t M 
for  0 £ t ^ T , ufe.L^jo,Tl  (that  is,  bounded  and 
measurable) , 

6.  the  condition  of  complete  controllability  of  chapter  1o. 
for  A,  F,  and  B be  fulfilled. 

Then  there  exist  a number  Ac  ^ 0,  a function  V(t)  of 
bounded  variation  on  [0,T]  and  a function  \j'(t),  solution 


- 4 - 


of 


. s 

■^(t)  = J(A*f  + Jf(x)F(c-s)dT  - Xe$x(xo,uo,s))ds  + V (t ) 


\<C 


on  [0,'f]  , such  that  X0  and  are  not  both  zero  and 
(B*'^(t)  - "Xo(l!IA.(xo,uo,t)  ,uQ(t)  ) = 


= max  (B%(t)  - *-0$u(x  .u.t)  ,u) 
utM  1 ° ° 


a.e.  on  |o,t)  . 
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3.  Some  historical  remarks. 


Problems  of  optimal  control  of  functional  differential 
equations  with  function  space  terminal  conditions  have  been 
investigated  only  for  the  last  years.  Among  the  approaches 
available  we  cite  the  followings: 

a)  BANKS  and  KENT  [8]  considered  the  equation 
■—  D(x ( • ) ,t)  = f (x(- ) ,u(t) ,t)  , 


D (x  ( • ) , t ) :=  x(t)  - J*ds  ;j.(t  ,s ) x (s ) . 

Vfi 

The  functions  are  taken  on  the  interval  [tQ-h,t^J.  This 
is  a neutral  functional  differential  equation  with  finite 
delay.  Using  Neustadt's  abstract  theory  of  extremals  sufficient 
and  necessary  conditions  for  minimizing 


v 

I = Jf°(x(t) ,u(t) ,t) it 


are  derived,  the  latter  in  form  of  an  integral  maximum  principle, 


b)  JACOBS  and  KAO  [29]  applied  the  abstract  Lagrange  multiplier 
rule  (see  LULNBERGER  [36]  ) to  nonlinear  delay-equations  with 
unconstrained  L^-controls. 

c)  BANKS  and  JACOBS  t6l  used  an  attainable  set  approach  in 
for  the  neutral  system 

x ( t ) = A.|  ( t ) x (t-h)  + A2(t)x(t)  + A3(t)x(t-h)  + B ( t ) u (t ) 

with  unconstrained  L2~controls  and  derived  a result  analogous 
to  those  in  a)  and  b) . 


d)  Recently,  BANKS  and  MANITIUS  [9]  used  the  projection  series 
method  for  linear  systems 
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x(t)  = L(xt)  + Du(t) 

with  L^-controls,  constrained  or  unconstrained,  and  quasi- 
convex  cost  functional.  They  obtained  the  optimal  solution 
as  the  limit  of  a sequence  of  solutions  to  certain  finite 
dimensional  problems. 

e)  DAS  [.18]  considered  a neutral  difference-differential 
equation  with  constraints  in  the  state  variable  and  obtained 
a generalization  of  b)  via  the  Dubovitskii-Milyutin  method. 

Applications  of  the  Dubovitskii-Milyutin  theory  to  general 
hereditary  systems,  such  as  Volterra  integrodif ferential 
equations  or  at  all  functional  differential  equations  which 
cannot  be  written  as  difference-differential  equations,  seem 
not  to  occur  in  the  literature  by  now.  The  reason  may  be 
complications  concerning  the  definition  of  an  "adjoint  system" 
in  connection  with  the  transformation  of  the  Euler  equation 
to  get  the  desired  maximum  condition. 

We  shall  use  the  space  C of  continuous  functions  as  phase 
space  and  L<^  -controls.  An  argument  against  L -controls  in 
the  theory  of  optimal  control  for  ordinary  differential  equations 
is  the  following:  If  the  range  M of  the  control  function  u 
consists  of  a finite  number  of  discrete  points,  it  is  difficult 
to  "variate"  the  control  function  in  a small  neighborhood  of 
uq  , due  to  the  Lt0  -norm  (see  GIRSANOV  [24]).  Now  in  view  of 
the  bang-bang-principle,  a finite  range  of  M seems  to  be 
"natural"  for  ordinary  differential  equations,  so  this  difficulty 
is  indeed  crucial.  But,  dealing  with  functional  differential 
equations,  this  argument  (at  least)  does  not  hold,  since  the 
bang-bang-principle  is  not  valid  (see  BANKS  and  KENT  [8]),  so 
a finite  range  of  the  control  variable  in  general  would  not 
allow  an  optimal  governing  of  the  system. 


4.  Definitions. 
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Let  S be  a given  function  space. 


S :=  \ 4>:  (a,b)  ->  RnJ  . 

I i j ' 


where  a or  b may  be  infinite,  and  x a function 

n 

X:  I — > R , 

I some  interval  in  R.  If  tel,  t-(b-a)  e I , then  we 
denote  by  xfc  an  element  of  S given  by 

Xfc  (9)  : = x ( t + 0 - b ) , © 0 ( a , b ) . 

If  t <L  I , t+(b-a)  6 I,  then  we  denote  by  xfc  an  element  of 

S given  by 

xt(0)  .=  x(t  + © - a) , © c (a ,b) . 


To  prove  the  main  theorem  we  need  the  definition  of  cone  and 
dual  cone.  Let  X be  a topological  vector  space  and  C <_  X. 

C is  a cone  with  apex  at  0,  if 

^ > 0 , x 6.  C »->  A x 6 C . 

Let  C be  a cone  v/ith  apex  at  0.  Denote  by  X'  the  dual  space 
of  X (the  dual  space  with  the  strong  topology  will  be  denoted 
by  X%<  ) . The  dual  cone  C*  is  defined  by 

C*  ;=  < f t X'  :Vx  & C:f  (x)  0 7;  . 


A nonzero  linear  functional  f is  called  a support  (or  supporting 
functional)  for  a set  A c X at  xq  A,  if 


x C A =>  f (X)  > f (xq) 
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^ • Auxiliary  mater ial . Adjoint  rclatior  s . 


We  shall  need  the  following 


Lemma.  Let 


L (xfc)  :=  j G (6) x(t+0)d9  , 


L*(^S)  :=  j'j’(s-e)  G (G) d© 


where  xq  = 0,  y = 0, 

variation  and  x continuous.  Then 
r T 

|'^(t)L(xt)dt  - j L*  (/)x(s)ds 


Proof.  Substituting  s = t + 0 we  get 

-r-  o 


•\tj  and  G of  bounded 


f 


- co 


j ^(t)L(xt)dt  = ^ y (t) ^ G (0)  x (t+0)  d0dt 

0 o 

T <■ 

('  f 

= jy(t)  G(s-t)x(s)ds 

0 — co 

T -r 


'iu(t)  G (s-t ) x (s)  dsdt 


iso  t 

*r  T 


= j Lj\i(t)  G (s-t)  dt . x (s)  ds 
i = o i.  ® !» 


fr! 


= , i vV  (s-0)  G (9)  d0  : x (s)  ds  = i 

rJ1  " J 

s O - «°  e 


1 L*(*ys)x(s)ds 


J 
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6 . Auxiliary  material . Existence  of  solutions  in  snace  V. 


Consider  the  equation 

T 

= j(A*^(s)  + L*'/)ds  + W(t)  , (6.1) 

i 


Here  the  n>;n  - Matrix  A';'(t)  is  continuous  and  of  bounded 
variation,  L+C''S  is  defined  as  in  chapter  5.,  and  W(t)  is  of 
bounded  variation  on  Lo,t]  . 

We  are  interested  in  solutions  of  (6.1)  in  the  soace  V of 
functions  of  bounded  variation.  The  following  is  valid: 

Theorem.  The  equation  (6.1)  has  an  unique  solution  in 
the  space  V[0,Tj  . 

We  prove  this  theorem  by  applying  the  contractive  mapping 
principle  to  the  operator 

T 

P v!>(t)  :=  (A*Us)  + L*V*)ds  + W(t)  . (6.2) 

i J » l 

i 

This  operator  maps  the  space  V[T-h,Tj  of  functions  of  bounded 
variation  with  0(T)  = 0 into  itself  for  small  h > 0 . The 
space  is  a Banach  space  with  the  norm 

r 

l!  9 II  :=  Var  n(t) 


We  use  the  following  relations: 


1 o 


lim  sup  (tk+l)  - A*(tk)|l  If  (tk+1)|J  + 


Him  sup{2lA’,(tk)Hlf(tk+1)  -lf<tk)(j 


Var  (A*  (t)  ) • sup^(t)  , T-h  4 t ■£  T j + 


+ sup^A* (t) , T-h  4 t 


T T 

4 £ 2M  Vai/lj>  (s); 


° T 

T t,  7 f T ( 

Var  (L  tp  ; = Var  ( (0  ( t-0)  G (©)  d©)  = Var  ( o (s)  G(t-s)ds) 

T-lv  1 T-t.  / ' T-f  J ' 

*'■  -to  1 w S.*t 


I I 

= lim  sup  2_||<f(s)G(tk+1-s))ds  - | ( j)(s)G(tk-s))dsJ  = 


= lim  sup  zlh 
- lim  sup^  J ( | tj> ( s )J  [ G (tk~s)j|  d 


-t»c 

T 


(s)G{tk-s)  )ds  + l/)(s)  (G(tk+1-s) -G(tk- 


fl f 

s + 


+ lim  sup2„  [ \({0 


z.1  R?(s)G(tk+1-s)  -G(t-s))dsj  . 
' i ' 


ktl 


Here  the  first  summand  is 


1 CP ( s ) J supjGjds  6:  h Vartp.  M , 

1 t-il 


T-L 


the  second  summand  is 


T 

4:  lim  sun/  j |?(s)!  |G(tk+1-s)  - G(tk-s){  ds 

T-(v 

T 

- lim  sud2. lG(tk+rs)  -G(tk-s)[)ds  ^ 

r-v. 


s)  ) ds 
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T 

^ h Var  (!9(s)  ) M. 

T-0.  1 


So  Var  (L*(pk  \ ^ 2hM. Var  ('j'(s)  ) 

1 T-f^  1 


3.  Let  f be  of  bounded  variation  and  set 

-fc 

F (t)  :=  J f (s) ds  . 

0 

Then 

•tk+i  “tf 


T 

Var (F (t) ) = 

0 


f f y 

= lim  sup  ^21.  | j f(s)ds  - j f (s)dsj^  = 


o 

ti«-i 

(V  I / _ v _ \ 1 


fc<ui 

r 


= lim  sup  j / | if  (s)ds|  ' ^ lim  sup  , Jf(s)[ds 

L J 

■t  t<  i-K 

T 

= |'|f(s)|ds  Tjjf(0)|  + Var(f(t))J  . 


By  the  same  consideration  we  also  derive 


T , T 

Var  (F  ( t)  ) ^ T,  f (T),  + Var(f  (t)  }.  . 

n *-  o ■* 


This  result  cannot  be  sharpened  as  one  sees  by  f(t)  = 1 on 

[0,1) , f (1)  = 0. 


Now  from  (6.2)  we  get 

T 

T r 

PV  - pyil  = Var(  (A*(-(s)  — ( s ) ) + L*(.i-^))ds 

T-*  ^ J 
t 


^ h Var  (A**(  (s)  - y(s)  ) + L*(jJ  - | ) ) 

*1- a 


« h(Var(A*(  (s)  - Y(s)  ) + Var(L'f(  -•/))) 

T-*.  v 1 11 


4 h(2M.  Var(  (s)  - Ws)  ) + 2hM  Var  ( (s)  -\f'(s))) 

TH.  1 1 T-fN  1 1 

= 2hM  ( 1 + h)  5 - Y !i  • 

Therefore  P is  a contractive  operator  for  h sufficiently 
small.  This  implies  the  existence  of  an  unique  solution  of 
equation  (6.1)  on  [T-h,T|  . 

Now  we  can  repeat  the  argument  on  [T-2h,T-h]  and  continue 
the  solution.  After  a finite  number  of  steps  we  have  proved 
the  theorem. 

Note  that  under  our  assumptions  there  exists  a constant  M 
such  that 

|A*(t)|  ^ M , Jg  (t)|  ^ M on  [0,T] 

r . 

Var  (A"  (t);  ^ M 

O 

T 

Var  (G (t) ) 4 M . 

o 


i 


r 


r\ 
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7.  Auxiliary  material . The  theorem  of  Dubovitskli-Mtlyutin 


The  theorem  of  Dubovitskii-Milyutin  provides  necessary 

O. 

conditions  for  the  solution  of^certain  extremum  problem.  This 
problem  can  be  formulated  very  general  so  to  enclose  a great 
variety  of  special  classical  and  modern  extremum  problems. 

To  formulate  the  theorem  we  introduce  the  following  notations. 

Let  X be  a locally  convex  topological  vector  space,  and 
xQfcX.  Consider  the  functional  F defined  in  a neighborhood  cf 
xq.  We  are  interested  in  a local  minimum  of  F(x),  where  x 
satisfies  two  types  of  constraints; 


(1)  x t ; i 1 , 2 , . . . , n 


where  C X,  4 !) 


(2)  x €.  C 


n+1 


where  Cn  + 1 <1  X,  Cn+1  = Cl  . 


Usually  , . . . ,C  are  called  inequality  constraints,  Cn+^  is 
called  equality  constraint. 

Now  let 

'h 

C ;=  Q C.  . 

We  loo!;  for  a x C such  that 

o 

F (x  ) = min  f (x)  , 

° C 

that  is,  we  want  to  determine  necessary  conditions  that  must 
hold  for  such  xQ. 

We  have  to  introduce  some  definitions. 

(a)  Directions  of  decrease . 

A h t H is  called  direction  of  decrease  of  the  functional  F(x) 


i 


at  xq,  if  there  exist  a neighborhood  U(h)  and  a c'.  = (!i) 

0,  such  that  for  0 S t < , h d.  U 

F(xq  + t h)  F ( xq ) + £(/. 

The  directions  of  decrease  generate  an  open  cone  K v/ith  apex 
at  0.  The  functional  F is  called  regularly  decreasing,  if 
the  set  of  its  directions  of  decrease  at  the  point  xq  is  convex 

(b)  Feasible  directions . 

A h & X is  called  feasible  direction  for  C.  at  x 

1 o 

(i  4 (l ,2, . . . ,nj ) , if  there  exists  a U (n)  such  that  for 
0 < t < t o , h t.  U 

x + e.  h 6.  C . . 
o 3 

The  feasible  directions  generate  an  open  cone  K with  apex 
at  0.  We  call  Cb  regular  at  xq  , if  the  cone  of  feasible 
directions  is  convex. 

( c ) Tangent  directions . 

A h C X is  called  a tangent  direction  to  C . . at  x , if 

n+1  o 

for  every  neighborhood  U (0)  there  exists  an  t0>  0,  that  for 

every  ' , 0 < t-  < f:  o , we  can  find  a x ( - ) c , such  that 

f r (£)  e U(0)  , 

C 

where  r(t)  :=  x(t)  - x - c h.  (In  a Banach  space  this  means 

!!  r (■• )(!  = ou  ) .) 


(Formulation  and  figure  in  GIRSANOV  p.  4o  could  be  misleading) 
The  tangent  directions  generate  a cone  K with  apex  at  0. 

Cn+^  is  called  regular  at  a point  xo,  if  the  cone  of  tangent 
directions  is  convex. 
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Now  the  main  theorem  can  be  given. 

Theorem.  Let  F(x)  be  regularly  decreasing  at  x^  and 

denote  by  the  cone  of  directions  of  decrease.  Let 

, i = 1,2,...,n  be  regular  at  xq  where  K.^  denote 

the  cones  of  feasible  directions.  Let  C , . be  regular  at 

n+1  J 

xq  and  denote  by  Kn+i  the  cone  of  tangent  directions. 

Let  F(x)  assume  a local  minimum  on  C at  x £ C. 

o 

Then  there  exist  K7  C.X'#  i = 0, 1 , . . . ,n+1  , not  all 

identically  zero,  such  that 


f 

o 


+ f. 


+ 


+ f 


n+1 


0 


(L) 


A proof  of  this  theorem  may  be  found  in  GIRSANOV  fi»d,p.41. 

In  applying  the  theorem  to  our  extremum  problem,  the  following 
tasks  must  be  solved: 

A.  Determination  of  directions  of  decrease,  feasible  directions 
and  tangent  directions. 

B.  Construction  of  the  dual  cones  K*  . 

1 

C.  Transformation  of  (E)  (the  "Euler  equation")  to  get  the 
desired  maximum  condition. 
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8.  Remarks  on  the  used  spaces. 


The  "true  state"  of  a hereditary  system  governed  by  a FDE 
is  an  element  of  some  function  space.  In  examining  control 
problems  with  FDEs  a principle  question  is  which  function  space 
is  to  be  used,  that  is,  how  we  have  to  topologize  the  set 
F = and  which  xfc  we  accept. 

In  contemporary  investigations  there  are  mainly  three  spaces 
that  have  proved  to  be  useful: 

a)  The  space  C of  continuous  functions, 

b)  The  Sobolev;  space  , 

c)  The  space  . 

Ad  a)  : As  to  the  domain  of  definition  we  ^iave 

C ;=  : trr  ,°i  ->  Rn  j with  the  sup-norm. 
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(See  BANKS  and  JACOBS  |X]  ) • 


Ad  c) : The  space  is  defined  as 

M2  : = Rn  XL2(  t-h,Oj  ,Rn)  , 


with  the  norm 


where  we  set  6=  ( <5>®,  6^  ) . (See  DELFOUR  ) 

( 1 ) ' 1 

The  Rt2  can  be  embedded  in  the  M2  in  a natural  way  .by  setting 


W2  :=  | ( <J»(0)  ,4)  fcM2#  e.  M2 


(see  MANITIUS  [£?_() 


In  applying  the  method  of  Dubovitskii-Milyutin  to  Volterra 
integro-dif ferential  equations,  two  points  have  to  be  watched. 

The  functions  x(t)  are  defined  on  (-o,T^.  That  is,  from  the 
preceding  spaces  only  the  C can  be  considered  as  an  appropriate 
space.  But  this  space  is  (only)  a metric  space,  not  a Banach 
space . 

One  crucial  point  in  the  application  of  the  method  is  the 
determination  of  the  cone  of  tangent  directions  to  the  equality 
constraint,  that  is  to  the  solutions  of  the  equation.  The 
determination  of  .this  cone  works  via  the  theorem  of  Lyusternik, 
which  assumes  Banach  spaces. 

Therefore  we  have  to  introduce  a new  space,  where  the 
following  considerations  take  place: 

1.  It  should  be  a Banach  space,  so  that  the  theorem  of  Lyusternik 
can  be  applied. 

2.  It  should  contain  the  continuous  functions  on  {-<*>, Tj  , since 
usually  system  states  vary  continuously. 

These  are  the  motivations  for  introducing  the  space 

C * (- «*,:'(:=  ^ h ••  (-  cO  ,r<  A Rn  , (|>  bounded  and 

continuous 

with  \\^>[|  = sup  - <*>  4 t & c^j  . 
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9.  Existence  and  uniqueness, 


In  this  chapter  we  show,  that  the  equation 


X ( t ) = A(t)x(t)  + j F(t-l)x(l)di  + B(t)u(t) 

(v)  J a.e.[o,T' 


Xo  = cj)0  c,  C^C-o/'I 

A ( t ) , B ( t)  continuous  on  j_0,Tj 

Y 

F(s)  of  bounded  variation,  jjF(s);ds  < co 


u 6 Lol°'Ti 

has  a unique  solution  on  [o,Tj . 

For  this  proof  we  use  a method  introduced  by  MORGENSTERN  \J\X y 
in  this  connection  first  used  by  STETTNER  [So],  The  method  is  based 
on  the  following  Lemma  by  CHU  and  DIAZ  : 

Lemma . Let  S be  an  arbitrary  set  and  A,K  mappings, 

A : S — > S , K : S S . 


Assume  there,  exists  a right  inverse  to  K,  that  is,  a Kr  , 

such  that  KK  ^ —I. 

r -! 

Then  A has  a fixed  point  in  S iff  AK  has  a fixed  point 
in  S . 


The  proof  of  this  lemma  is  almost  self-evident. 

To  prove  the  unique  existence  of  a solution  of  (V)  in  the 
space  of  bounded  continuous  functions  on  (-<ro,'fj  , we  introduce 

the  operators  A and  by 


<Mt) 


- s£>  t £ 0 


Ax  ( t ) ; = j -fc  ;v 

^ 4>,(0)  + j (Ax(s)  + j F(s-t)x(t)dt  +Bu(s))ds  , 

P - Co 

O £ £ £ T 
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( x(t)  - CO  < t £ O 

K,x(t)  : = 

A I It 

I e x(t)  0 i t ^ T . 


-1 


Kk  has  an  inverse  operator,  namely  KA 


K^x(t)  = 


/ 

_ J 


x(t)  - co<  t £ O 


e“  Atx(t)  0 £ t £ T. 


,-V 


Now  consider  the  operator  AK^ . For  two  elements  x,y 

c'3(-  '•i  ,Tj  we  get 

Hk_1AK  x - K_1AK  y ||  ^ 


£ sup  |e  ?>t  |(AKax  - AK^y)  (t)  | ; 0 ^ t ^ T ] 


, l 


$ 


sup  |e-/U:  j [e’'s  M^jx  - y !/  + 1} x - ylfv  j e'  C 


o 

O i t 


Here 


We  used 


t 

je"XsM1i!x  - y|J  £ A - y|(eXt  , 
o • 


•fc  s 


o o 


U(t)\  L Mr  |f(s)[  ^ V . 

Now  the  sup  is 

(.  sup  ^m1  A”1 1!*  - yll  + v-x"2!^  " y ft  j 

& M2  X” 1 U x - y|| 


of 


dtjds  , 


1o.  The  cone  of  tangent  directions. 


Consider  the  equation 


x(t)  = Ax ( t ) + Lxfc  + Bu(t), 


Uo-\) 


where 


Lxfc  ;=  f F ( t-  ~c ) x ( - ) dr.  , 


with  the  boundary  conditions 


co  = e - 


r = § C,  C^T-r ,Tj ^ 


(T-r  > 0) , 


under  the  hypotheses  of  chapter  9.  This  equation  defines  a certain 
set  C . = v(x,u)  in  the  space  C*"  x Lto  . We  want  to  characterize 

the  cone  of  tangent  directions  to  the  set  C , . 

J n+1 

Introduce  the  operator 

P:  C X Lw  *->  C **  x C [ T-t,  T] 


j'  x (t)  - j>0  (t)  , -c»  < t £ 0 


P(x,u)  : = 


I ! • 

y x ( t ) - (■.•„(())  - j (ttx  (s) +Lxs+Bu  (s)  ) ds  , 


0 <tiT  ; 


' XT  ^ 


Then  the  solutions  of  equation  ( 1 o . 1 ) are  the  zeros  of  the 
operator  P,  that  is,  they  satisfy 


P(x,u)  = 0. 
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To  find  the  cone  of  targent  directions  to  the  solutions  of  the 
last  equation,  we  use  the  following  theorem  of  Lysternik  (see 
LYU  STERN  IK  and  S0130LEW  Q3b])  *. 
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Theorem . Let  X,Y  be  Banach  spaces,  P an  Operator 


P:  X -»  Y 


vanishing  for  xy  G X.  Assume  the  Frechet-derivative  P'(x) 
exists  and  is  continuous  in  a neighborhood  of  xq.  Let 
P'  (x  ) be  surjective. 

Denote  by  K the  cone  of  tangent  directions  in  to  the 
set  Q = x ; P ( x ) = Oj  . 

Then 


K = £ h : P ' (xQ)h  = 0 j 


To  apply  the  theorem  of  Lyusternik  to  our  problem  we  first  show 
that  the  operator  P has  a continuous  Frechet-derivative. 


We  see 


P(x+x,u+u)  - P ( x , u ) = 


x (t)  , 


- *©  / t £ 0 


• x ( t ) - . Ax  ( s ) + Lx  + Bu ( s ) ) d s , 0 M t £ T 

\ J s J 

O 


Since  the  right  side  id  a linear  continuous  operator,  we  get 


P' ( x , u ) (x,u)  = 


x (t)  , 


- <o  < t k 0 


) - / (s)  + L>:  + Bu  ( s ) ) ds , 0 


; x,r 


■t  - .•  r i ••at  i vt  , 


namely  the  validity  of 

P(xQ  + h)  - P(xo)  = P 1 (xq) h + r(xQ,h), 

v/here  |l  r [!  = o(|!hl\).  In  our  case  r ~ 0 . Note,  that  P is  not 

linear ! 

The  continuity  of  P'  follows  by  known  arguments  (see 
LUENBKRGER  , p.193). 

Now  we  want  to  show  that  P'(x,u)  is  surjective.  That  is, 
for  arbitrary 

((a(t),bT)  €.  C*°(-«o#T'i  x C ['f-r/i'j 

we  must  solve 


x(t)  = a (t) 


- r t £ 0 


(S)  x(t)  = a(t)  + (^Ax(s)  + Lxg  + Bu(s))ds  0 £ t ,£  T 

© 


This  system  has  a solution,  if  the  condition  of  complete 
controllability  is  valid.  Therefore  we  assume 

(CCC ) : The  System  (S)*  is  completely  controllable. 

Then  we  can  apply  the  theorem  of  Lyusternik  and  get  the  following 
result: 

The  cone  of  tangent  directions  to  the  set  is  given  by 


K:=|X,u:x  =0,  x(t)=  i (Ax(s)  + Lx  + Bu(s))ds, 

V O y s 

o 

O £-  t £ T,  xT  = 0*  e.  C « ,Ta  X La,[0,Tj 


The  condition  of  complete  controllability  means  the  following: 
Consider  the  system  (S) 1 
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■ y(t)  = Ay  (t)  + Ly.  + bu(t)  O ^ t T 

I 

I 

j 

(s)  • •!  y(t)  = o,  - <o  a ^ o 

u & Lo\0,Tj 

As  shown  in  chapter  9.,  this  system  has  a solution  on  [o , T j . 
Denote  by  F the  set  of  all  yri  6 C [T-r , Tj  , where  y(t)  is  a 
solution  of  (S)  ' , as  u ranges  over  all  of  Ij0oCc>/Tl* 

(S) ' is  called  completely  controllable  if 

F = c[T-r,Tj  . 

Now,  if  (S) 1 is  completely  controllable,  (S)  has  a solution.  To 
prove  this,  chase  z(t)  such  that 


z(t) 
z (t) 


a(t)  , - co  < t £ 0 


(t)  + j(. Az(s)  + Lz^ds 


0 £ t £ T , 


(this  has  a solution  by  chapter  9.),  and  then  chose  y(t)  such  that 


y(t)  = 0 , -p'j  < t £ 0 

y(t)  = j(Ay(s)  + Lys  + bu(s))ds  0 < t^  T 


i 

I 


(this  has  a solution  by  the  condition  of  complete  controllability) . 
Now, 


x (t)  :=  y (t)  + z (t) 

u ( t ) ; = u ( t ) 


is  a solution  of  (S)  . 
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11.  The  dual  cone. 


We  want  to  determine  the  dual  cone  to  the  cone  of  tangent 
directions  derived  in  chapter  1o.  Since  K is  a subspace  of 

X :=  c‘°  (-^Tj  y LtojO,T^  , 

we  get  (see  GIRSANOV  { , p . 6 9 ) 

K ^ =^ffcX,r:  f (x , u)  =0  for  (X/iileKj  . 

Let  us  characterize  the  functionals  in  K*  . To  do  so  we  use 
a lemma  about  adjoint  operators  (see  LULNBERGLR  £2(4,p. 1 56) : 

Lemma . Let  X,Y  be  Banach  spaces  and  T a linear  operator 
T:  X— > Y,  bounded.  Assume  that  the  range  of  T be  closed 
Denote  with  T*  the  adjoint  operator  of  T,  with  ...  (T*) 
range  of  this  operator,  with  _/.r(T ) the  nullspace  of  T 
with  S'1  , SCX,  the  "orthogonal  complement",  that  is 

Sl  :=  ^x*  £ Xv‘ : x*  (x)  = 0 for  all  x £ Sj 

Then 

fc(T*)  =i^f(T):X  . 

*■  *i 


Now  define  an  operator 


: c ~ (-cvrj  * L J0,T)  ->  Y C 'jl’-r  ,T 


by 


T(x,u)  := 

■ x ( t ) 

I 

^x(t) 


- ® < t i 0 
t 

(Ax(s)  + Lx  + Bu(s))ds,  0 

J s 

o 


t 


(actually  this  is  the  operator  P'(x,u)  derived  in  chapter  1 
and  denote  elements  of  C (-  to,T,  x L ^0,T  shortly  by  x 


the 

and 


J 


I 
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Then 

K*  * Lf  <L  X*:  f(x)  =0  for  Tx  = 0 ’ , 

that  is,  x is  "orthogonal"  to  f,  if  x //(T)  . 

Now  7C(T)  is  closed,  since  the  system  is  completely  controllable. 
Therefore  we  can  apply  the  previous  theorem  with  the  following 
result: 

f must  be  in  the  range  of  T * , that  is,  there  exists  a y, 

* f6* 

y € [c*5*  (-»,Tj  x-  c[T-r,Tjj  . = : D 
such  that  f = T*y.  Hence 

f (x ) = (T*y) (x)  = y(Tx)  . 

Now  the  elements  of  the  dual  space  D have  the  form  (d,0) , 
where  0 is  a function  of  bounded  variation  on  [T-r,T]  (to  be 
more  precise,  the  second  coordinates  form  a space  isomorphic 
to  the  space  of  functions  of  bounded  variation) . The  form  of 
the  first  coordinate  need  not  be  known,  as  soon  will  be  seen. 

If  we  choose 

x (t)  = O - *t  S 0, 

o 

•fc 

x(t)  = i Ax ( s ) + Lx  + Bu(s))ds  , 0 ^ t £ T , 

■J v s 

o 

then  the  first  summand  in  y(Tx)  vanishes  and  we  get 
f(x)  = <x,r,0^>  , 

where  , •/  denotes  the  "pairing"  of  elements  of  C^T-r,!’  with 
elements  of  its  dual  space,  that  is, 

T 

f (x)  = j x(s)d0(s)  , 

T-r 

the  latter  being  a Stielt jes-integral . 


L 


i 


PH"".'  | V-  - ' 

1 2 . The  Euler  Liquation  . 


Now  consider  our  control  problem  (CP) . Let  x ,u  be  a solution 

o o 

°f  this  problem.  We  shall  establish  the  Euler  equation  by  applying 
the  theorem  of  Dubovitskii-Milyutin.  To  do  this,  we  have  to  determine 
the  cone  of  decrease  and  his  dual  cone  for  the  functional 

T 

F ( x , u j = j (£>(x(t)  ,u(t)  ,t)  dt  . 

0 

Using  a known  result  concerning  functionals  of  this  type  (see 
for  instance  GIRSANOV  , p.56,  Example  7.7  and  p.69,  theorem 

1o.2)  we  get; 

If  K f O , then  for  any  f £ K * , 

o'  Joo 


T 

fo(x,u)  = “ j*‘(i>x(xo,uo,t)  ,x)  + (<!>u(xo,uo't)'u))dt' 


where  > 0. 


Next  we  analyize  the  constraint 

u(t)  t M for  0 £ t £ T a.e . 

Denote  by  the  cone  of  feasible  directions  on  the  set  at 
xo'uo'  where 

Q1  = C x Q1  , 

Qj  c the  set  of  functions  u(t)  with  values  in  M a.e. 
Repeating  an  argument  of  Girsanov  m , p. 86 , b) , we  get ; 

If  f^  4.K*  , then  f^  = (0,f  * ) , where  f.|  <L  is  a support 

to  Q.J  at  uQ. 

Now  we  apply  the  theorem  of  Dubovitskii-Milyutin  to  (CP)  and 
derive: 


There  exist  functionals 
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f .f,,f,  ^ 

O 1 2 


not  all  zero,  such  that 


fQ(x,u)  + (x,u)  + f2(x,u)  = O , 


(12.1) 


where  f ,t.  just  were  defined,  and  f9£K*  (see  chapter  11.),  that 

O 1 


f 2 (x)  = x (s)  d©  (s) 


xQ(t)  = 0 


-fO<t  ^ O 


:)  = j(Ax(s)  + Lx's  + Bu(s))ds,  O 2 t ^ T 


(12.2) 


Then  (12.1)  becomes 


f 1 (u)  = Aoj 


T 

(xQ,uo,t),x)  + f<[\  (xo,uoft)  ,u)  )dt  - 


- I 


x (s) de (s) 


(12.3) 


In  the  following  we  transform  the  right  side  of  (12.3)  by  using 
the  function  V defined  by 

T 

\}/(t)  = [(A*^+  A<,<f>y)ds  + e(T)  - e(t). 


O £ t £ T 


= 0,  (especially  i|»(T)  = o)  . 


(12.4) 


Here  we  assume  e(T-r)  = 0 (that  is,  0 is  normalized  with  respect 
to  the  Stielt jes-integral  in  (12.3),  see  PFLAUMANN-liNGKR  ) 


0(t)  = 0 for  O £ t < T - r. 


In  chapter  6.  wo  proved  the  existence  of  solution  to  this  equation. 
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Now  we  first  want  to  show  that 

.T  T J 

AoJ  ( 4>x  *)dt  - jXde  = - j (13*^,u)ds 

° T-f  ° 

The  right  side  of  this  equation  is 

T 


- j(B*  ^ ,u)dt  = - j (V^,Bu)dt 


o 

T 


= - j ( , x - Ax  - Lxt)dt 


(12.5) 


(12.6) 


We  shall  shov/ 

T 


j((^,X  - Ax  - Lxt)  + ( '*)  ) dt  = 


T--r 


T 

= jxd©  “ Y^T-r  ^ x (T-r^ 
T-'f 


T-f 


f 

( (^,x  - Ax  - Lxfc)  + (\.'^,x))dt  = 


0,  T-T 


(12.7) 


= ^(T-r)x(T-r) 


(12.8) 


Addition  of  (1^.7)  and  (12.8)  using  (12.6)  proves  (12.5) 
Nov/  with  (12.4)  we  get 
T 

j(.(  ,x  - Ax  - Lxt)  + ( )l#r^,x))dt  = 


TW 


T T 


= j(  j (a*^  + £ \ - ),4^)ds  + e(T)  -(e(t)),5c)  - 


T-C  t 


- Ax  - ‘^Lxt  + ;uix)dt  = 

r T 


= Jxj(A*^  + L^-^ldsdt  - }6Xdt  + 0(T)  ^ xdt 

T-r  * T“* 


- j ( ''■['Ax  - \;Lxt  + A</^x  x)dt  = 


T-* 


3o 


I 


= Cj(A>  + L y - ACJ<v)  ds,x  (t)  j 


T-i 


f 


- j-(Air^  + L*  f - A0^,x)dt  - 

T-r 

+ e(T)x(T)  - e(T)x(T-r)  - 


r 

j ©xdt  + 
T-f 


f 


J 

T-f 


( Y Ax  - Y Lxfc  + A6's>sx)dt  = 


l 

= - ( | (A*  y + L* ^ - A.<$>x  ) dt , x (T-r ) ) - 


T-f 


+ 0(T)x(T)  - e (T)  x (T-r)  = 


\ 0xdt  + 

j 

T-r 


-(  (T-r)  - 0 ( T ) + 6(T-r)  ,x  (T-r) ) - 0x|  + 

T-t 

T 

+ jxd©  + © (T ) x (T ) - 0 (T)  x (T-r)  = 


\ 
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k I 


ft  I 

I 


i I 


r 


j ( ^ Ax  + yLxfc  - lo^xjdt  = 


T 

( | (A*^>  + hK^'-%0fyy  )ds,x(t)  ) 
b 

T-r 

c 


T-< 


T-f 


- j - (A*  y + L*^-  A.  <^x  ,x)dt  + j 6 (T)'idt  - 


T-f 

r 


- j ( y,Ax)  + yLxt  - Xo^x^)  dt  = 


(since  x(0)  = 0) 


( j (A*4'  + L*^-  >-.^x  ) dt ,x (T-r ) ) + ©(T)x(T-r) 
T-f 


= ( y(T-r)  - 6(T)  + ©(T-r)  ,x  (T-r)  ) + 0(T)x(T-r)  = 


= y(T-r ) x (T-r) 


This  proves  (12.5).  Nov;  (12.3)  becomes 

T 

f'  (u)  = f (~B%(t)  + >;6(x  .u  ,t)  ,u(t)  )dt  (12.8) 

1 J > iu.  O O 


-/ 


with  arbitrary  u , where  f ^ (u)  is  a support  to  at  uo> 

We  use  now  the  following  lemma  (see  Girsanov  [3*1  J , p.7G,  example 
1o.5) : 


Lemma.  If  the  linear  functional 

T 


a (t ) 6 L1 (0,T) 


(u)  = \ (a  (t)  ,u  (t)  ) dt  , 

o 

is  a support  to  Q*  at  uQ  , then  u£  M implies 


(a  (t)  ,u  - uq)  >•  O 


i.e.  [0,t]  . 


With  this  lemma,  we  get  from  (12.8); 
u fe.  M implies 


I 


13.  Remarks. 


The  preceding  result  is  not  unexpected.  Indeed  it  has  the 
usual  form  of  such  maximum  principles.  It  is  similar  to  the 
results  cited  in  chapter  3.  or  to  the  "classical"  results  for 
ordinary  differential  equations  (see  LEE  and  MARKUS  (341). 

Though  maximum  principles  of  this  kind  are  satisfactory  from 
a theoretical  point  of  view,  for  the  explicite  computation  of 
the  optimal  control  function  only  half  work  is  done.  The  next 
step  must  be  the  numerical  evaluation  of  the  condition.  This 
evaluation  is  possible  for  special  cases  of  ordinary  differential 
equations  and  shall  be  investigated  for  Volterra  integro- 
differential  equations  in  the  next  part  of  this  work. 
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